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ABSTRACT 


This thesis investigates the problem of automatic 
transcription of the morse signal, and describes and 
documents several approaches to filtering, processing, and 
decoding it for transcription. The baseband signal is first 
modeled as a modified random telegraph wave. A discrete 
Kalman filter and a linear smoother are then used to process 
the demodulated signal in order to gain a measure of the 
effectiveness and applicability of this model. It is shown 
experimentally that this model and processing yield a 
Significant reduction in the transcription error rate. Next, 
a Viterbi decoder algorithm based on a Simple Markov model 
of the code is programmed and tested. Finally, the base- 
band signal model is incorporated in a more general model 
for pre-detection Kalman filtering. It is shown that this 
filter permits acceptable recovery of morse signals whose 
average signal-to-noise ratio is as low as -14 dB in a 2 kHz 


bandwidth. 
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т. ЖЫРКОРОСТТОЧ 


Economic inflation and the national commitment to the 
all-volunteer Armed Forces concept have combined to produce 
unprecedented increases in costs of both manpower and weapon 
systems in recent years. The public's keen awareness of 
these higher costs, together with less than enthusiastic 
support of defense programs in general, has caused Congress 
to be reluctant to authorize increases in defense expenditures. 
Thus it has become necessary to reduce the number of armed 
forces personnel in order to keep defense expenditures within 
authorized limits. This reduction has had the effect of 
intensifying the development of mechanization of appropriate 
manual tasks on a broad front. 

Signal surveillance conducted by the armed forces, 
recognized as an essential and integral part of intelligent 
tactical and strategic planning, is one such area where 
automation is receiving increased attention and support. In 
particular, the human operator has long been relied upon to 
provide the necessary manual transcription of manual morse 
circuits under surveillance. Because of the reduced manpower 
levels, this surveillance and transcription must be tranferred 
to mechanized equipment if this source of intelligence is to 
remain timely and effective. 

This thesis investigates the problem of automatic trans- 


cription of the morse signal, and describes and documents 
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several approaches to filtering, processing, and decoding 

it for transcription. The baseband morse signal process is 
first modeled as a modified random telegraph wave. A discrete 
Kalman filter and a linear smoother are then used to process 
the demodulated signal in order to gain a measure of the 
effectiveness and applicability of this model. It is shown 
experimentally that this model and processing yield a 
Significant reduction in the transcription error rate. Next, 
a Viterbi decoder algorithm based on a simple Markov model 

of the code is programmed and tested. Finally, a more general 
model of the signal process, incorporating the baseband model, 


is used to design and implement a pre-detection Kalman filter. 
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ТТ. PROBLEM DESCRIPTION 


A. THE MANUAL MORSE SIGNAL PROCESS 

It is assumed that the reader is familiar with the 
manual morse signal, its pecularities, vagaries, and 
uncertainties, and with current methods of transcription. 
To formalize the discussion, however, certain definitions 
of the terms used and a brief description of the signal are 
in order. 

Пе ОИ пошто ENIS report, the term morse signal 
Be eho кета оспа! Моге Ссае, sent manually by key, 
manual "bug", or electronic keyer. The problem of tran- 
Sem bingekey board morse, that which is sent automatically 
with standard parameters, will not be considered, although 
рано тиште Се аге арріпсар1е. The baseband morse signal 
is the output of the keyer and 1s represented by the logic 
levels "0" and "1", corresponding to the states "key up" 
and "key down." The five characters of the international 
morse code are identified as: dot, dash, element-space, 
letter-space, and word-space, The term element re to 
the standard time unit of the code; its actual duratıon in 
seconds will of course vary with sending speed. Standard 
morse code consists of the character lengths shown in 


Table I, 


eS tne terms baud ana bit are used. 
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TABLE I 


STANDARD MORSE CHARACTERS 


Character Symbol Duration (in elements) 
Dot 5 1 
Dash 3 
Element-space А 1 
Letter-space E 3 
Word-space и 7 


The standard word (including word-space) in morse 
communication is 50 elements in length. Thus the element 
duration in seconds for a given sending speed may be 
calculated as 6/5 times the reciprocal of the speed in words 
per minute. The author is unaware of any generally accepted 
standard for the bandwidth of the baseband morse signal; it 
was found to be convenient to express the upper limit of the 
bandwidth as three times the reciprocal of the element 
duration. Thus a code speed of 36 wpm has a bandwidth of 
SU HZ. 

An actual (as opposed to standard) morse signal, as 
those familiar with the problem are aware, may exhibit 
quite a wide variation from standard code in character 
duration, speed variability, and consistency of element 
duration. Since these variations are often unique to the 


particular sending operator, and in many cases may depend 
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ӘРНЕГЕ ог Еғгатгіс being sent as well, it is difficult, 
if not impossible, to describe a "typical" morse signal. 

This variability is illustrated in Figures 1 and 2, which 
are histograms of the character duration of two different 
amateur ratio operators recorded on the air. As can be seen, 
the distributions are different, with the least variability 
in both cases appearing in the dot and element-space 


durations. 


DU SYSTEM CONSTRAINTS 
The signal processor will obviously play an important 
role in the automated manual morse intercept and transcription 
system. In order to fully appreciate the system constraints 
under which the processor was developed, and the context in 
which it 1s expected to operate, an outline of the integrated 
system is presented. Referring to the system block diagram, 
Figure 3, its three basic components may be briefly described 
as follows. 
1. Modulation Sorting Subsystem 
This subsystem scans the band and/or frequencies of 
interest and detects the presence of morse signals in a 
(typically) 2 kHz band. Upon detection of a morse signal, a 
Separate digitally-tuned receiver is automatically tuned to 
the signal frequency for reception. 
2. Morse Processor 
This signal processor is currently the only part 


of the system which is not in existing hardware or software. 
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орао function 15 to minimize error probability in the 
face of noise, interference, and uncertain signal parameters. 
3. Automatic Transcriber 
This component translates the code into letters of 
the alphabet; there are currently several transcribers 
available which have proven effective at adequate signal- 


to-noise ratios with modest signal parameter variation (11,121. 
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MIT Про РО OBJECTIVES 


The ‘ultimate goal of the intercept and transcription 
system is to provide output copy with an error rate no 
greater than that which a "good" human operator can provide. 
Such an operator manually transcribing the morse Signal can 
adapt rapidly to changing signal parameters and has little 
trouble distinguishing dots from dashes even if the sender's 
keying is far from perfect. Additionally, he can adapt 
readily to the noise and interference environment and reliably 
copy a Signal in the presence of numerous other morse and 
non-morse signals. 

Specific operator performance data were not available to 
the author, thus only broad design objectives were formulated 
based on a limited number of both subjective and experimental 
data obtained using amateur radio operators as subjects. 
Random letter sequences were sent using the Pickering model 
КВ-1 погзе keyboard to key a signal generator at an audio 
frequency selected by the subject. Noise was added to the 
audio signal and the SNR in each bandwidth used was recorded. 
The results, summarized in Table II, tabulate error rate 
versus SNR in the bandwidth used. Also shown is the SNR in 
the signal bandwidth as previously defined in Section II.A. 

The conclusions drawn from these data is that a good 
operator can copy reasonably well down to -13 dB SNR in a 
2 kHz bandwidth. Although the insertion of a 100 Hz bandpass 


filter raises the SNR to O dB, the relative invariance of 
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TABLE ІІ 


OPERATOR PERFORMANCE DATA 


(a) Speed: 35 wpm 


BW SNR SNR ERROR RATE (%) 
(ав) (ӘВ) 
(in given EW) (in signal BW) ОРІ OP2 
2 KHZ = 3 0 12 15 
OO HZ | => 0 11 13 
100 Hz 0 0 12 14 
ZEKHz -10 3 6 10 
200 Н2 0 3 11 9 
100 Hz 3 3 11 10 
2 kHz -7 6 1 2 
200 Hz 3 6 0 1 
100 Н2 6 6 2 2 
(b) Speed: 25 wpm 
BW SNR SNR ERROR RATE (3) 
(dB) (dB) 
(in given BW) (in signal BW) ОРІ OP2 
Рана, -13 2 6 5 
200 Hz -3 2 7 6 
100 Hz 0 2 8 5 
2 kHz --10 5 2 1 
200 Hz 0 5 1. 1 
100 Н2 3 5 2 2 
2 HZ E 8 1 0 
200 Н2 3 8 0 1 
100 Hz 6 8 0 2 
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difficuit and 
fatiguing 
BW too narrow 
fatiguing 
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relatively easy 
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COMMENTS 


difficult and 
fatiguing 
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Dut still milaly 
fatiguing 


easy enough 





error rates with changes in filter bandwidth indicates that 

the ear performs the necessary filtering. Strict concentration 
is required, however, at this low SNR, and the test operators 
stated they would not attempt to copy such a signal unless 
strongly motivated. 

Using the previously defined signal bandwidth for each 
speed, the results may be summarized as follows: An operator 
can provide copy with a 10-15% error rate with an SNR of 
approximately 0 dB in the signal bandwidth; at 3 dB, the 
error rate is 5-10%; and at 6 dB, copy is practically perfect. 

Based on these results, it seems reasonable to assert 
that a typical field operator, faced with searching for and 
copying morse traffic eight hours a day, would not be inclined 
to copy signals much below 6 ав SNR in the signal bandwidth 
unless absolutely required. Thus the system designer must 
select from two costly alternatives: 1) If he designs the 
system to perform as well as the good operator is able to 
perform, the automated system will reliably receive a large 
percentage of signals encountered on the air, but it is likely 
to be complex and expensive. 2) On the other hand, if he 
Шесте Еле wsystem tO perform as well as a typical operator 
probably performs, the automated system will be cheaper, but 
the remaining operators who must copy the low SNR signals 
which are not machine transcribable may become too fatigued 
to be efficient, leaving the overall man/machine surveillance 
system less effective than the existing manual system. 

Such design considerations are beyond the scope of this 


thesis; using the "good" operator as a criterion, the 
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ultimate system design objectives may be broadly stated as 


follows: 


1) 


2) 


With an error rate of 10% or less, recover and 
decode morse signals whose SNR in a 2 kHz bandwidth 
is on the order of -10 dB, using standard code with 
additive white gausSian noise and no interference. 
Track the time-varying statistics of character 
lengths in order to enable the transcriber to 


translate the code properly. 
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N ROSES SUE DESIGN PHILOSOPHY 


The first step in the processor design was to model the 
morse code as a random telegraph wave with non-stationary 
transition probabilities. Using this model, several 
increasingly complex processing methods were implemented, 
and the processing gain of each stage was determined. First 
a Kalman filter was designed to filter the demodulated 
output of a square-law detector with and without narrowband 
analog IF filtering. Next a smoothing algorithm was added 
to determine the effectiveness over Kalman filtering alone. 
Finally a maximum a posteriori (MAP) estimator of the code 
characters, using the smoothed output for likelihood 
calculation, and using the Viterbi algorithm for processing, 
was programmed and tested. 

After determination of the error-reduction effectiveness 
of each of these processing stages, a more general model 
of the signal process was used to design a Kalman filter for 
pre-detection filtering. The objective was to determine 
whether or not such filtering yielded any advantage over the 
Simpler demodulation/post-detection filter approach. A block 
diagram of the various stages is shown in Figure 4. Sections 
V through VII present a theoretical basis for each of the 
processing stages, followed by a presentation and discussion 


of experimental results. 
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V. BASEBAND MODEL AND PROCESSOR 


The baseband morse signal may be modeled as 
x(k+1) = x(k) + w(k) (V-]) 
where x = key state (0 or 1), and w(k) is a random forcing 


function descriptive of the morse keying process. At the 


output of the demodulator, the Signal is observed as 
Z(K) = X(k) +t (К) 


where z(k) is the observed value and v(k) is the additive 
noise. This model gives rise to the following scalar 


Kalman filter algorithm [3],[4]: 


V(k|k-1) + R 
V(k|k) = [1- G(k) ]V(k|k-1) (estimation variance) 
у(к+1 |к) = V(k|k) + Q(k) (prediction variance) 


x(k|k) = х(к|к-1) + с(к) Га(к) - x(k|k-1)] 


(estimation) 


х(к+1 |к) = х(К|К) (prediction) 
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where: 


x(k) = estimate of x at time k 
R = variance of the observation noise, v(k) 
Q(k) = variance of the random forcing function, w(k).. 


Since this algorithm implies knowledge of the variances R 
and Q, procedures for their estimation are required. In 
order to keep the filter algorithm TT as simple as 
possible, estimates of Q and R were made independently, and 
used by the filter algorithm as if these were the true values. 
A. ALGORITHM FOR ESTIMATING THE RANDOM 

FORCING FUNCTION VARIANCE 

The random forcing function w(k) is descriptive of the 
on-off keying process, i.e., it may be thought of as the 
process which causes the transitions in the x state from 
0 to 1 апа from 1 to 0. Referring to the signal model 


equation (V-1), the keying process has the following 


interpretation: If x(k) = 0 and w(k) = 0, then x(k+1) = O 
and x remains in the space condition. If, on the other hand, 
w(k) = 1, x shifts (at time k+1) from the space condition to 


the mark condition. This process is illustrated in Figure 5. 
Since the probability of occurrence of a transition is 
dependent on the time duration since the last transition, 
this probability is non-stationary. A proper description of 
the transition probabilities at each time k, then, must 


necessarily be time dependent and conditioned on the element 
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Figure 5. 
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Illustration of Signal Model Process 
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duration апа the present state. The following probabilities 


are therefore required: 


Р (м (К) =0 |К,х(К) =1) Pr[remain in mark condition] (V-2) 


P(w(k)=1|k,x(k)=1) (У- 3) 


||! 
O 


P(w(k) =-1[k,x(k) =1) Pr [transition from mark to space]  (V-4) 


P(w(k)=0|k,x(k)=0) = Pr[remain in space condition] (V-5) 
| P(w(k)=1|k,x(k) =0) = Pr[transition from space to mark] (V-6) 
P(w(k)=-1|k,x(k)=0) = 0) (V-7) 


Probabilities (V-3) and (V-7) are identically zero since 

state values less than 0 and greater than 1 are not allowable. 
The remaining probabilities are dependent on the distributions 
of dot, dash, element-space, letter-space, and word-space 
durations of the particular morse signal being received, and 
are dependent in a Markov sense on the previous character. 

The Markov nature of the code character transitions was not 
taken advantage of in the filtering process. Additionally 


it was assumed that a particular operator's character 


27 





durations are all distributed with uniform density, that 
the dot and element space have the same density, and that 
the dash and letter space have the same density. The 
assumed densities, as shown in Figure 6, presume for the 
present that the mean values T, Та, and the parameters %, а 
are either known or have been determined in some manner. 
Although the assumption of uniform densities for 
the character durations of any particular sending operator 
is probably not strictly justifiable, neither is the 
assumption of any other well-known density, such as a 
gausSian or exponential density. The most likely candidate 
for properly modeling these duration distributions may be 
a gaussian-like density with the tails truncated at suitable 
values. The complexity of estimating the parameters of 
such a density for a particular received signal, together 
with the computational burden of evaluating the error 
Accion) ЕГІ(ЕС), for the probability calculations (V-2) 
through (V-7) at each sample point, motivated the selection 
of the uniform density. The resulting probability computations 


are relatively simple and straightforward. 
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a) Dot/Element-space Duration Density 


Рас Ê) 
E. 
2%, 
0 t 
Т-2 T T+2 
b) Dash/Letter-space Duration Density 
Pdash‘®) 
] 
2d 
: t 
142% 1a T ytd 
Figure 6. Character Duration Densities 





l. Mark Transition Probabilities 
A description of probabilities (V-2) and (V-4) may 
be obtained by first conditioning on the dot and dash 
probabilities, and noting that if x(k) is still in the mark 
condition after T+, then the mark is known to be a dash. 
Given then that a mark, m, is a dot, and given T and ў, 


probability (V-2) may be modeled as 


T+ $ 
À » ET 
Py (dot) = P(w(k)=0|k,T,x(k)=1,m=dot) = : A CHE 
lo; o к <= 
24 T-k | B 
= эў + 5 і E See Se 
0; Шак 
Similarly, for a dash: 
тата 
Py, (dash) = P(w(k)=0|k,T¿,x(k)=1,m=dash) = : Paasn(t) dt 
i Ш 57,5 0% тата 
Thk 
= а | _ 
= as Та ТЕТЕ тата 
о; Та+а < к 
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Here k is the time index which counts the number 
the signal has been in the mark condition on the 
mark. The conditioning on m must be removed for 


since m is not known during this interval. Thus 


Of Е5 
Present 
ОСЕК < ШИРА 


for k < Т%%: 


P(w(k)=0|k,T,x(k)=1) = Py (dot) Pr(dot) + Pq; (dash) Pr (dash) 


апа Гог К > Т+% 
P(w(k)-0[k,Tg,x(k)-1) = Po; (dash) 
But the dot and dash probabilities are dependent 


of traffic in the message, 1.e. on what language 


is in, whether it is plain text or code groups, 


(У-8) 


(У-9) 


оп the type 


the message 


letters only 


or both letters and numbers, etc. Since the traffic type 


may not be known a priori, equiprobable dots and 
assumed at this point. Using Pr(dot) = Pr(dash) 


equations (V-8) and (V-9) reduce to 


1 У ss 

D + i Пе < 

P(w(k)=0|k,T,T,,x(k)=1) = 1 T+ < 
Hh + š аша 

0 Tp td < 
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dashes were 


= 1⁄2, then, 


~ 
IA 


Т+ 


~ 
IA 
+3 





Probability (V-6) follows immediately as 


IO ТАО ЕН) Е (w(k)=0|k,T,T,-(k)=1) 


The expression for probability (У-2) ав а function of time, 
k, conditioned on the present state and dot and dash 
distributions is sketched in Figure 7. 

псе i Gansitlon Probabilities 

An appropriate description for the space duration 

probabilities (V-4) and (V-7) is derived similarly, first by 
conditioning on a particular space and then removing the 
appropriate conditioning by using the relative frequencies 
of each space. Given that a particular space, s, is an 


element-space, then, 


1 en S 
x à T- l 
P(w(j)=0|x(j)=0,s=elem,T) = ЭЙТУ T 
0 TFL < j 
(V-10) 


where j is the time index which counts the number of units 
the signal has been in the space condition on the present 
space; and the element length, given T, is assumed to have 


the same uniform density as the dot length. 
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Similarly for a letter-space: 


P(w(j)=0|x(k)=0,s=ltr-sp,T,) 


1 0 <j < TQ, 

т 

= 5 T, Soy таға (Vall) 
0 таға I M 


where j is the same index as above and the letter-space 
length, given Та! has the same density as the dash-length. 

For the word-space it was decided to use an 
exponential model, since after about 5T units, the word-space 
is (in actual practice) about equally likely to end at any 


time: 


P(w(3)=0|x(3)=0,s=word-sp,T) = 





(v-12) 
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The conditioning on s may be removed by observing that if 
je таға, the space must be a word-space; if 3 > Т“2, then 
the space is a word-space or letter-space. In equiprobable 


letter traffic, the needed probabilities are: 


Pr(s-elem-space) - 12/17 
Pr(s=letter-space) = 4/17 
Pr (s=word-space) - 1/17 


Applying these probabilities then to (V-10), (У-11), апа 


(V-12) yields the desired expression: 


P(w(j)-0|j,T,T4,x(3)-0) = 





1 0<34<т-4 
p) Z2 Э. M 
1 THE < j « T--d 
1 =) 
а аа -а< < т.+а 
ББ 14452244 
1 Таға IS 
– (4257) 
2T 9T « j 
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Dx catculatrion of Variance 
The variance of w(k) is now easily calculated at 


cael LINC, Ku 4S follows: 


Me I Mr "ee En дил ut 


Let 


р, (0,1) £ Pr(w(k)=0|x(k)=1) 
ру (-1,1) £ Pr(w(k)=-1|x(k)=12) 


Pr(w(j)=0|x(j)=0) 


|| 


un 


Pr(w(j)=1|x(j)=0) 


|| 


Р. (1,0) 


Hen, k)w-» 1, 


2 


(0)*-P, (0,1) + (-1)*+P, (-1,1) - [E(w(k)) 1]? 


0 (к) 


2 
EX) = са ери 


Py (-1,1)-P, (0,1) 
Similarly, if x(j) = 0, 
Ola meee (0,0)-P:(1,0). 
J J 
Since, at the receiver, the true state x(k) is not 
known, correct estimates of Q are dependent on correct 


estimates of the state, x(k). Thus, it is expected that at 


some SNR, incorrect estimates of x(k) will cause the 
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estimates of Q to be erroneous enough to force incorrect 

estimates of x(k+1), causing a runaway condition to develop 
and yield the receiver worthless. This SNR at which runaway 
develops was determined experimentally, and found not to be 


a serious problenm. 


ЕРІ СПАБАСТЕБК DISTRIBUTION ESTIMATION 

It was assumed that the character distributions are 
appropriately described by uniform esito: with known 
parameters (expect for mean values). More sophisticated 
methods of distribution and/or parameter estimation were 
discarded in favor of simplicity. Experience suggests that 
even asloppy sender will usually not exceed an 2/Т or d/T 5 
ratio of about 1/3. Thus, once estimates of T and Та аге 
obtained, £ and d can be determined from this assumption 
unless these ratios are known in advance. The mean values 


Ш апа Т of dot or element-space and dash or letter-space, 


а 
respectively, were estimated by measuring the character 
durations and sequentially averaging the appropriate duration. 
Thus T is the mean value of the dot and element-space durations, 
and Та is the mean value of the dash and letter-space durations. 
TNE EnNresNOld fOr deciding which set of measurements a 
particular length belongs to was set at the halfway-point 
between dot and dash lengths. The algorithm is as follows: 
Initially specify: 


t; = shortest dot-duration expected, 


СЕ 
|| 


longest dot-duration expected, 
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1. (а) Initialize estimate of T as (t4*t5)/2. 
(b) Initialize estimate of Та as (3t,+3t,)/2. 

2: -Measure"the'duration-of "each mark and space. 

3. If the measured value is less than (3t,+t,)/2, then 
identify it as a dot/element duration, Ti. 

4. If the measured value is larger than (3t,+t,)/2, then 
identify it as a dash/letter-space duration, T,. 


5. Estimate the means recursively by 


1 


(ai T (K) = T(k=1) + i [Ti = ү к=] 


(b) Talk) = Tq(k-1) + E (7, - T, (k-1)] 
С. OBSERVATION NOISE VARIANCE 

Although the noise power used in each of the tests was 
known and could have been specified initially, it was 
decided to estimate this parameter in order to better 
Simulate an operational environment where the noise power 
is not known a priori. Although the method is purely 
intuitive and without a valid theoretical basis, reasonably 
good results were obtained. 

At the output of the square-law demodulator, the noise 
1s no longer gaussian and is correlated with the signal. 
Proceeding, however, as if the signal and noise were not 
correlated, the noise variance R can be obtained by sub- 
tracting the morse Signal power from the total received 


(demodulated) Signal power: 
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~ 


R A 22 
ағыу саты” 


where Me is obtained recursively by the following: 


У, (К) = V,(k-1) + 2 [(2-u,)% - V,(k-1)] 


mean value of the received (demodulated) signal 


o> Ea 
|| 


estimate of demodulated signal amplitude. 


~ 


The parameters түй and a were also estimated on line, although 
they too are known a priori for test purposes. The estimator 


algorithm for а 15 


[2(К) - а, (к-1)] if zik) > u, 


жік 


1) а. (К) а. (К-1) + 


[z(k) - aj(k-1)] if 2(K) < Wy 


^im 


(2) a,(k) 2 a4 (k-1) + 


(3) а(К) = а, (к) - а, (к). 


The mean value, ШТ is simply 


u. (k) = u,(k-1) + È [z(k) - u,(k-1)] . 
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D. SMOOTHING ALGORITHM 

The Kalman-filtered output represents the optimum 
linearly-filtered estimate of the signal amplitude in the 
minimum mean-square-error sense based on the assumed model 
and the efficiency and consistency of the estimated 
statistics. This estimate may be improved, however, by 
smoothing the data, which implies consideration of future 
inputs as well as past inputs [4]. 

The equations for the discrete optimal linear fixed- 
interval smoother, for the scalar case, expressed as a 
combination of a forward-running filter and backward-running 


filter are given Бу [5],[6] аз: 
Smoothed Estimate: 
x(k|N) = P(k|N).[x(k|K)/Pz(kK|K) + x (k|k+1)/P (k |k+1)] 


Estimation Variance: 


1 


MU - IZP-UETE)- Ж І?Р, (ЕТЕРІ) | 


where 
x(k | k) = filtered state estimate 
P.(k|k) = error variance for the forward filter 
x, (k|k+1) = (predicted) estimate for the backward filter 
Py (kK[k*1) = (predicted) error variance, backward filter. 
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These expressions are not computationally suitable, however, 
since P, (k|k+1) is not finite when k = N. Alternate 
expressions [6] which are more amenable to computation are 


given as follows: 


P (kk) 


E и (К) = т A ау "° 


Треп 


P(k|N) = (1 -W(%)/P, (k]x+1)]%-P,¿(k|k) - [6% (x) /P, (k]k+1) ] 
and 


1/P, (k|k+1) = 0 when К = КМ 


Let W, (k|[xk+1) = x, (k|k+1) /P, (k|x+1) 


Then 


а | х(к |к) 
х(к м) = тт + P¿(x]x) /P, (K]k+1)7 + P(k]N)w, (X[k+1) 


and 


P(k|N)W, (k|k+1) = 0 when k =N 
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These equations then represent the smoothing algorithm. 
As noted above, forward-filtered estimates must be stored 
until the same data can be backward-filtered and combined 
to produce the smoothed state estimate. 

The smoothed state estimate then is the best (linear) 
estimate possible for the given model and parameter estimators, 
and any improvement in error rate must be derived from giving 


more probabilistic structure to the signal model. 


Е. IMPLEMENTATION OF FILTER AND SMOOTHER 

A diagram of the test Signal generation method is shown 
in Figure 8. The square-law demodulator was selected simply 
for ease of implementation since an analog Squarer and 
appropriate filters were readily available and easily inter- 
faced with the analog-to-digital converter. The 100 Hz 
low-pass filter permits recovery of morse signals of 
approximately 35 wpm or less. The signal-to-noise ratios 
used throughout this report are average (pre-detection) 
Signal-to-noise power ratios and not pulse signal power to 
noise power ratios. More specifically, the average signal 
power in the morse signal is defined as 

P. = >= Т [s(t) ]^ at 

where the interval [0,2T] is a dot and element-space. This 


expression reduces to 
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Р_ = sa / (a cos wt)? dt + > Í Оша 
0 EN 


P_ = а2/4 ; 


which is 3 dB less than the pulse signal power. 


The noise signal is taken from a calibrated white gaussian 
noise source of 1 volt rms. 

The filter and smoother, along with the auxiliary 
estimation ECL ms) were coded in Fortran and implemented 
on the XDS-9300 computer interfaced with the CI-5000 analog 
computer for analog filtering and for D/A and A/D conversion. 
The sampling rate was 500 samples per second, and the value 
of N for smoothing was chosen to be 250 samples. The sampled 
test signal was recorded on tape for subsequent processing, 
since the processing required approximately 4 seconds for 
l second of data. The test signal runs consisted of the 


following: 


1. Perfect code AR sequence and random letter sequence at 
Speeds of 35, 30, and 25 wpm each with a signal-to-noise 
ве 5 4, 3, 2, and 1 dB in a 2 kHz BW with no 
pre-detection analog filtering. 

2. Perfect code AR sequence and random letter sequence at 
speeds of 35, 30, and 25 wpm each with a signal-to-noise 
ratio of -7, -8, -9, -10, -11, -12 dB in a 2 kHz BW with 


a 100 Hz pre-detection analog bandpass filter. 
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3. Sloppy code AR sequence at nominal speeds of 25 and 30 
wpm each with a signal-to-noise ratio of 6, 5, 4, 3, 2, 
1 dB in a 2 kHz BW with no pre-detection filtering. 

4. Sloppy code AR sequence at nominal speeds of 25 and 30 
wpm each with a signal-to-noise ratio of -7, -8, -9, 
-10, -11, -12 ав іп а 2 kHz BW with a 100 Hz pre-detection 


analog bandpass filter. 


The output of the processor was recorded on an 8-channel 
strip-chart, using a utility D/A conversion routine. Тһе 
outputs were as follows: 

Channel 1: Recorded input signal 


Channel 2: Unprocessed output (input signal thresholded 
at its mean value) 


Channel 3: Kalman filtered output 

Channel 4: Filtered output thresholded at 0 V. 

Channel 5: Smoothed output 

Channel 6: Smoothed output thresholded at 0 У. 

Channel 7: Option of filter gain or noise variance estimate 


Channel 8: Dot/element-space duration estimate. 


The following figures (9-19) are typical output records, 
showing examples of the test runs for each signal and type 
of sequence. Channels 1-4 are shown in Figure (a) in each 
case with the corresponding channels 5 and 6 shown in Figure 
(5). An example of the output of channel 7 for the gain 
option is shown in Figure 20, along with the corresponding 


outputs from channels 1, 3, and 4. Figure 21 shows an 
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example of the outputs for the noise variance estimate апа 

the dot/element-space duration estimate at a point where the 
input signal changes from 25 to 30 wpm and the signal-to- 

noise ratio changes from 1 dB to 6 dB. In all cases the 

chart speed was 10 mm/sec except in Figure 21 where 5 mm/sec 

was used in order to show the estimates more clearly. The scales 
for channels 1, 2, 4 and 6, and for channel 7 gain option, 

are 5 v./div., with 100 volts corresponding to a variable 

value of 1.0; the scales for channels 3 and 5 are 2 UAIN: 

For the channel 7 variance option the scale is 200 mv./div., 


and for channel 8 the scale was calibrated at 4 msec/div. 


BE RESULTS ӨК TESTS 

The outputs of the processor were decoded by hand to 
determine the error rates. Using the estimate of T, a mark 
was decoded as a dot if its duration was 2T or less and as 
a dash otherwise. Similarly, a space was decoded as an 
element-space if its duration was 2T or less, as a letter- 
space if the duration was between 2T and 4T, and as a word 
Space otherwise. 

The following data were obtained for each run: 
1. Letter error rate and bit error rate with no processing. 
2. Same error rates with filtering only. 
3. Same error rates with filtering and smoothing. 

A letter error occurs when any transmitted letter is not 
correctly decoded. Only one error per transmitted letter 


is counted; for example, if "A" is decoded as "ET", one 
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letter error has occurred. Оп the other hand, if "ET" is 
decoded as "A", two errors have occurred, since neither "E" 

nor "T" was decoded correctly. Word spaces are counted as 

one letter per 7 elements. The sample size in each case was 
approximately 200 letters. A bit error is defined as at 

least one mark-space error occurring within a transmitted 
element duration. Again, only one error per element is counted, 
and the sample size was approximately 200 bits. 

The results of this analysis are presented in Tables III 
through VI. Column 1 lists the signal-to-noise ratio used 
without a 100 Hz pre-detection filter; column 2 lists the 
Signal-to-noise ratio used with the 100 Hz filter in place. 
Both the bit and letter error rates for filtering and 
smoothing are tabulated, with the error rates for the 
unprocessed output shown for comparison. In each case, the 
error rates are shown for the typical (random letter) 
sequence and the AR sequence, except in Table VI, where the 
results are for the AR sequence only since no random letter 
sequence for this case was recorded. Table VII shows 
typical hand-translated sequences for each processing stage. 

These results indicate that the Kalman filter and linear 
smoother provide a Significant decrease in both bit and 
letter error rates. By using a 100 Hz bandpass pre-detection 
filter, such processing provides a tolerable 10% letter error 


rate оп а -7 ав SNR signal as opposed to an unacceptable 
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Smoothed Output for 35 wpm AR Sequence, 4 dB SNR 
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Smoothed Output for 35 wpm Typical Sequence, 4 dB SNR 


FIGURE 130. 
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TABLE VII 


TICE TRANSLATED SEQUENCES 


SENT 
UNPROCESSED 
FILTERED 
SMOOTHED 


SENT 
UNPROCESSED 
FILTERED | 
SMOOTHED 


SENT 
UNPROCESSED 
FILTERED 
SMOOTHED 


SENT 
UNPROCESSED 
FILTERED 
SMOOTHED 


SUMMARY : 


total letters sent 
unprocessed errors 


Mon) 0 KO 
e) e) n Je. 
ل۰ ل ل‎ 


QAZ 
IEQAZ 
IEQAZ 
EEQAZ 


QAZ 
QUAD 
OMZ 
QAZ 


QAZ 
OC EZ 


QAZ E 


ОА2 


filtered errors 
smoothed errors 


(question mark indicates 


BW = 2 kHz 


SNR = 5 dB 


WS ХЕ D 


WH X E 


VV VF 3 
nun 
[т] 
дж» 
HHH 


Wore E E 
асан ы I 
WHXEEI I 
WHX II 
WS X E 
LS X EIH 
WS X EEE 
WS X Е 


€ 


ЕЕЕ ЕТЕР? 
Е ИН X IEDC 
рс 


DE 
DC 
DY 
DC 


De 
2» 
ІШЕ 
р? 
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KE VT EB YHN 
E TEV TEG E THN 
Е ВРУ T GG YHN 
IV ОВ YHN 


De RFV T 
DCERE ER 
ШЕ НАУ 
ЕТ 


GB YHN UJM 
СВЕУНМЕО?? 
GB YHN UJM 
GB YHN UJM 
REV T GB YHN UJM 
Bese ve LT 
REV T 
REV T 


GB 
GB QUNEUJM 
YHN U 
Е ҮН?ЕМ 
SINT U 
MENU 


Ju 
LFVIA 
REVET 
КРУ Т 


СВ 
СВ 
СВ 
СВ 


Error Rate 
46% 
25% 
15% 


untranslatable sequence) 


speed = 30 wpm perfect code 


UJM 
FIM 
UJM 
UJM 


PBEYHNE?JME 
UNES E 





32-34% error rate without processing. Smoothing contributes 
considerably to a reduced error rate down to about -9 dB 
where its effectiveness over filtering alone begins to fall 
off. For the high speed 35 wpm signal, runaway of Q estimation 
occurs at approximately -9 dB, while for the lower speeds 
runaway never really develops until the SNR reaches -12 dB, 
since filtering always provides an improvement in error rate 
even at these low SNR's. 

It was noted that a majority of the errors in the filtered 
and smoothed output result from insertions of isolated dots 
in the letter-space and particularly the word-space separations. 
A possible remedy to this situation is to incorporate the 
Markov structure of the code in the estimation algorithm for 
О, although an increased susceptibility to runaway may limit 
its effectiveness. 

Since the bit error rate must be on the order of 1% or 
less in order to yield a tolerable 10% letter error rate, 
it was felt that bit error rates of up to approximately 3% 
could possibly be reduced to an acceptable level by use of 
soft-decision Viterbi decoding following the smoothed output. 
If such a reduction were possible, then acceptable error 
rates could be obtained for SNR's down to -9 dB. 

Additionally, pre-detection Kalman filtering and parameter 
estimation would yield a theoretical gain of 3 dB or more 
over non-coherent demodulation. With these improvements in 


the processor, then, the output letter error rate of 10$ 


0 





could be achieved on signals whose SNR is on the order of 
-12 dB in a 2 kHz bandwidth and the originally stated design 
objective would be met. 

Both of these possible improvements were implemented 
separately to determine the effectiveness of each one. The 
Viterbi decoder algorithm used the smoothed output of the 
post-detection Kalman filter as input. A separate program 


was written to implement a pre-detection Kalman filter. 
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NU VITERBI DECODER 


The Viterbi algorithm [7], [8], as originally formulated, 
is a maximum likelihood (ML) sequence estimation decoding 
algorithm for convolutional codes. It has found application 
in other areas [9], however, and its use has been extended 
to maximum a posteriori (MAP) estimation. It is the MAP 


estimation use which is of importance here. 


A. MAP ESTIMATION 

In order to transform the smoothed output of the 
processor into characters of the morse code, certain 
decision criteria must be formulated. The easiest and most 
obvious way to accomplish decoding is to threshold the 
smoothed output at its mean value and determine the identity 
of marks and spaces on the basis of the measured duration 
of each received character. Such a scheme, however, fails 
to utilize two sources of information which are inherent 
in the smoothed output: 1) Thresholding discards all 
information present in the actual smoothed amplitude estimate, 
and 2) the Markov nature of the character transitions 15$ not 
utilized. 

The decoding problem, then, is to take advantage of this 
additional information to determine the most probable 
sequence of morse characters. The thresholded output may be 
выса tO Make tentative decisions to obtain a specific 


sequence of character outputs ن2‎ 11 where a particular 


19 





2. 15 either a dot, dash, element-space, or letter-space. 
It is the objective, then, to find a specific sequence of 
transmitted characters Хе X, which maximizes the 


probability that X r+... X, Was transmitted, given that 


Zoro +. 12, Was received. Formally, it is desired to find 
the sequence х; = х. Хр which gives 
Р(хех.|2=а.) = maximum 


where ЕЕ 15 the specific received sequence, and all possible 
nin c so transmitted characters, x, are used to determine 
which specific sequence yields the maximum value. The 
sequence = which yields the maximum is then the maximum 
likelihood estimate of the transmitted sequence if it is 
assumed that the transmitted characters are equiprobable; 
it is the maximum a posteriori (Or minimum error) estimate 
if the actual.probability of transmission of each character 
шэ utilized [10]. 

In general, it would be necessary to compute and compare 
mae probability, Р(х=хр|2=2.), for all possible sequences x,. 
However, if it is assumed that the morse code is a Markov 
source, then the problem of finding the MAP estimate reduces 
to a problem of maximizing a sum and the Viterbi algorithm 
may be used. 

In the following development, a shorthand probability 
notation is used to facilitate writing of probability 


statements. The statement P(x,|z,) is used to mean 
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P(x, =a; [2,=a5) where the a,,a, are the characters of the 
code, i.e., dot, dash, element-space, letter-space. No 
confusion should result, since in all cases the notation 
P(u) is intended to mean the probability that u is equal 


to some specific value. 


В. SOURCE MODEL 

A third-order Markov model of the morse code exhibits 
a good deal more probabilistic structure than a first-order 
model, as can be seen by comparing the transition probabilities 
shown in Table VIII. In the interest of simplicity, however, 
it was decided to use a first-order model. The assumption of 


a first order model then means that 
P (Xp | Xor ee e r Xk) E P (x [x 1) 


where Xy is the kth character of a transmitted sequence. 


The transition matrix lists the following transition 


probabilities: 
P (xy, [x i) for each x, = a, 
where 
a, =. (dot) (VI-la) 
Eg ol (dash) (VI-1b) 
ај = ~ (element-space) Ке) 
a, = ~ (letter-space) (УІ-14) 
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TABLE VIII 


MARKOV TRANSITION PROBABILITIES 


ENEST-ORDER MODEL 
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(Word spaces were not considered as separate characters, 


but as combinations of letter-spaces and element-spaces.) 


E SEQUENCE PROBABILITIES 
Based on the first order model assumption, the probability 


of any particular transmitted sequence of length n is given by 


ne 
Р(х ,ху,...,х _у) = I P (x; |x;.1) (VI-2) 
i-0 
Then, given an input sequence to the decoder (Zx °° 2 1), 


use of Bayes' rule expresses the desired conditional 


probability as 
Р (хои. _1|2,...,2,_,) 


E Ое Ze ах 3) P(X ore ees X 
DONE 7 ) 
O n-1l 


which is the probability to be maximized. 
Assuming that the thresholded output is memoryless, the 


conditional output sequence probability becomes: 
P(Z ге ер] | Хо XL) == I P(z; |x;) н (УІ-4) 
Although the thresholded output is by no means memoryless, 


due to the decision directed nature of the O estimation 


algorithm in the preceding filter, this assumption 15 
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nevertheless made in order to render the computation of the 
necessary probabilities tractable. The dependence of the 
thresholded character output decision on previous and future 
decisions will be removed to a certain extent in Section D, 
justifying the assumption of memorylessness at this point. 
Then, using (VI-2) and (VI-4), and realizing that 
Р(2,...,2 _1) Еее ан Е, maximization of (VI-3) is 
equivalent to maximizing the expression: 
1 Du 


P (x; |х; 1) Tl Р(г. |х.) i 
0 1= 


п 


ПЕРІ! 


3 


Maximization of this sequence is equivalent to minimizing 

the negative logarithm, since 11 Р(-) is a monotonic function 

ee P(-). Thus 
L(X r+. X 1) - i Па P(x,|x; 4) + In P(z, |x,)] 

is the function to be minimized by the Viterbi algorithm. 


An outline of how the Viterbi algorithm performs this 


minimization is presented in Appendix B. 


р. LIKELIHOOD COMPUTATION 
The likelihoods P(z,|x,) may be computed from the smoothed 
Signal amplitude and received character duration as follows. 
Define the following figures of merit for amplitude and 


duration: 
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ТЕ 


values 


спин се сасе of merit for a dot-length duration 


time figure of merit for a dot-length duration 


amplitude figure of merit for a dash-length duration 


time figure of merit for a dash-length duration 


these figures of merit are scaled such that their 


are between 0 and 1, then they may be interpreted 


es follows: 


an > 


Probability that 
interval. 


probability that 
probability that 


Probability that 
interval. 


probability that 


робар су that 


а 


а 


а 


mark occurred during a dot-length 


mark is a dot, or 
Space is an element-space. 


mark occurred during a dash-length 


mark is a dash, or 


space is a letter-space. 


Likelihoods, then, may be computed by utilizing these values 


as probabilities. 


For example, the probability that the 


thresholded output is a dot, given that the input 15 ап 


element-space is simply: 


Р(2=а. | х=а-) 


ta (AT) a, 


Pr(z=dot|z=mark) -Pr(z#dash |z=mark) -Pr(z=mark) 





The likelihoods for each case were determined as above 


and are given as follows: 


P (z=a, |x=a,) Zee) 


P(z=a,|x=a,) = a (1 - Є (1-7 )) 


P(z=a,|x=a,) EL uu S 


|| 


P(z=a,|x=a,) BE (0 - t (er )) 


for 1 - 1,3 where the а. are given by (VI-1); and 


P(z—a, |x=a,) Е ut) 


|| 


Р(а=а. |х=а.) пет =) 


EM тае) 


шас Sd) (1 2)1 (1€) 


Bor i = 2,4. 


The figures of merit a, A, t , T were obtained by 
m m m m 
using the thresholded, smoothed output to make "tentative" 
decisions, and then computing the merit of these decisions. 


The tentative decisions were: 
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ШІ ТІ ве = and measured duration < 2T, then x, = dot 
2) If x = ] апа measured duration > 2T, then x = dash 
out 

3) If xj, = 0 and measured duration < 2T, then x, = element- 
space 

4) If en 0 and measured duration ? 2T, then = letter- 
space 

For x, = dot or element-space, the previous T values of 


smoothed output were averaged and the amplitude figure of 


merit was taken as 


a eos 
0 = 2 > х. (Т) 
а x А S xc Be 
а + x, (T)/a 9O05 5 2 > > 
-- а 
1 ШЕ 


where 
х. (Т) = the average of the smoothed output, x,, over the 
previous T seconds. 
a = amplitude of smoothed output signal. 


The time figures of merit were 


1 0 < Са = T 
CA 2- (Е +1)/2Т T « t4 « 3T 
0 т = ta 


87 





0 0 < +4 < 74/3 
mom ЕЕ (morta) ТЗ ЗЕ; < Та 
1 | Ta < ta 


where ta - measured duration of character duration. 


These functions are sketched in Figure 22. For s = dash 
or letter-space, the previous Та values of x, were averaged 


to obtain the value for An: 


NIY > 


D 
|| 
м 9» 


х. (Та) /а + 0,5 


NW > 


iL х (1) > 


These likelihood computations allow decisions to be made 
by the Viterbi algorithm to determine the most likely 
character on a character-versus-character basis. It was 
necessary, however, to extend these computations to cover 
more complicated situations such as that depicted in Figure 23. 
The thresholded output shown decodes as .4.4. , although the 
sequence .4— was actually transmitted. The Viterbi algorithm 
implemented using the above likelihood computations decoded 
the sequence as .ı-ı. , however, since no provision vas made 
to account for the non-memoryless nature of the thresholded 


output. 
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a) Amplitude Figure of Merit 
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b) Duration Figures of Merit 
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Figure 22. Sketch of Figures of Merit 
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Smoothed 


Output 
0 y \/ Time 
- | 
Thresholded 
Output 
0 RUTA === au Tas t i -. 
Eine 


Figure 23. Illustration of Signal Requiring 
Modified Likelihood Computation 
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Two modifications were made, therefore, in order to 
account for situations such as this: First, when a dash or 
letter-space was decoded by the thresholded output, this 
character was divided into three equal segments and dot and 
element~-space likelihoods were computed for each segment. 
Secondly, when the thresholded output consecutively decoded 
three short characters, then the total duration of all three 
characters was obtained and dash and letter-space likelihoods 


were computed for the total interval. 


Е. IMPLEMENTATION AND RESULTS 

The algorithm was coded in Fortran and run as a 
subroutine of the original Kalman filter and smoothing 
programs. Sequences of 9-11 elements in duration were 
decoded. In order to determine the contribution of context 
information supplied by the first order transition 
probabilities, the algorithm was also run with equiprobable 
transition probabilities, i.e. as a ML sequence estimator. 

The results are summarized in Tables IX and X. Again 
both bit and letter error rates are presented with the hand- 
decoded smoothed output repeated here for ease of comparison. 
Table XI shows a comparison of bit error rates for ML and 
MAP estimation showing a slight improvement provided by the 
MAP estimation, taking advantage of transition probabilities. 

The reduction in error rate, although not as great as 
hoped for, shows that the thresholded output bit error rate 


can be improved significantly at the higher SNR's, but 
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TABLE ІХ 


VETERBI DECODER ERROR RATES - 35 wpm 


AR SEQUENCE, PERFECT CODE: 


SNR (2 kHz) SMOOTHED OUTPUT VITERBI OUTPUT IMPROVEMENT RATIO 
IN BIT ERROR RATE 
NO PRE-FILTER BIT LTR BIT LTR 
(cb) (9) (3) (9) (9) 
б 0.83 9 „35 3 2.4 
5 0 27 2.5 23 158 
4 6.7 49 6.8 49 = 
3 35 96 35 100 = 


mye ECan SEQUENCE, PERFECT CODE: 


SNR (2 kHz) SMOOTHED OUTPUT VITERBI OUTPUT IMPROVEMENT RATIO 
IN BIT ERROR RATE 

NO PRE-FILTER BIT LTR BIT DIR 

(db) (%) (%) (%) (%) 

6 0.70 8.0 0.33 3 DL 

5 1.2 15 0.90 Dp 103 

4 4.5 34 4.4 34 = 

3 30 92 30 96 = 
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TABLE X 


VITERBI DECODER ERROR RATES - 25 wpm 


AR SEQUENCE, PERFECT CODE: 


SNR (2kHz) SMOOTHED OUTPUT VITERBI OUTPUT IMPROVEMENT RATIO 
IN BIT ERROR RATE 

NO PRE-FILTER BIT LTR BIT LTR 

(db) (%) (%) (%) (%) 

6 227 10 ШЕСІ 2 25 

5 p 17 10 pos 

4 1 7 18 125 17 pon 

3 4.5 36 4.4 38 = 


ПЕРСА SEQUENCE, PERFECT CODE: 


SNR (2kHz) SMOOTHED OUTPUT VITERBI OUTPUT IMPROVEMENT RATIO 
IN BIT ERROR RATE 
NO PRE-FILTER BIT LTR BIT LTR 
(db) (%) (%) (9) (5) 
6 0.70 8 и 3 19 
5 D 12 0.80 10 1.4 
4 eZ 15 0.90 13 13 
3 29 24 2] 22 m 
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TABLE Al 


COMPARISON OF ML AND MAP ESTIMATES 


SNR (2kHz) MAXIMUM LIKELIHOOD MAXIMUM A-POSTERIORI 
MO PRE-FILTER ESTIMATOR ESTIMATOR 

(db) BIT ERFOR RATE BIT ERROR RATE 

(5) (%) 

6 0.45 0,5? 

5 Doom 0.80 

4 HE 0.90 

3 2.8 207 


Typical sequence 25 wpm, perfect code 
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152 
1 
2 
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ЕШ егког rates ot about 3% or higher remain practically 
unaffected. Again, it was noted that most errors are 
isolated dot insertions during word-spaces. Since the third 
order model accounts for the small transition probability 

of such events, such a model should virtually eliminate this 
type of error, resulting in a very significant decrease in 
error rates. 

Table XII shows a typical decoded sequence as output by 
the algorithm, with the filtered and smoothed outputs shown 
for comparison. This sequence was part of the 35 wpm AR 
sequence at 5 dB SNR in 2 kHz. Since it was noted that the 
algorithm created some new errors as well as correcting 
errors, likelihood computations could probably be improved 
by employing better character-length density estimation 
procedures. 

Two examples of likelihood computations and the evolution 
of the most likely sequence as the algorithm progresses are 
shown in Tables XIIIa and b. Table XIIIa is an example of 
decoding with highly probable likelihoods, while Table XIIIb 
shows an example of a correction by the algorithm. The array 
of numbers headed by SURVIVOR SEQUENCES indicate to which 
previous node the node at stage k is connected; the length 
of each survivor sequence is shown immediately below this 
array. The line through the survivor sequence array shows 
the final path for the minimum-length sequence. The next 


array shows the computed log-likelihoods for each character, 
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followed by the minimum length sequence with its length. 
(The number 200 was sufficiently large to be used for 
infinity.) Finally, the mark and space durations as 
determined by the thresholded output are shown for comparison. 
In these tables (XIIIa and b), the following numbers 


correspond to the decoded morse characters: 


1 - dot 


tN 
| 


element-space 
3,4 - letter-space 


5,6 - dash 


The letter-space and dash are represnted by two numbers since 
such decisions may result either from comparisons made on a 
straightforward character versus character basis or from 
comparisons made utilizing the modified likelihood computation 
previously described. In each case the sequence А. А—А.~ 

was transmitted; Table XIIIb shows how the Viterbi algorithm 
corrected the thresholded sequence (4.4.4.^4.-.) to form a 

dash from the inner ... sequence. This situation is similar 


to that shown previously in Figure 23. 
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EU DTEÉERED 


SMOOTHED 
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COMPARISON OF VITERBI OUTPUT WITH 


SMOOTHED AND FILTERED OUTPUTS 
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1. Sequence of AR at 35 wpm, 5 db SNR, tranmitted. 


SUMMARY OF VITERBI OUTPUT: 


New errors made 

Errors uncorrected 
Corrections 
Net improvement 
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S ORE DETECTION MODEL AND FILTER 


Since it is assumed that the receiver has been tuned to 
the signal carrier frequency, it 1s sufficient to model the 
Signal as one of known frequency and unknown amplitude and 
phase Of the form A sin(wt+0). The received signal then is 


of the form 
AS A БИ (ЕТО) T cv (t) 


where v(t) represents additive white gaussian noise. A state 
model of the signal process may be obtained by rewriting the 
transmitted signal in the form a sin wt + b cos wt [4]. 
Letting ху represent the discretized Signal in this form, the 


following state model results: 


x, (k+1) jk WT ху (К) 

хо (К+1) ZWT il x, (k) 
where X, (t) B b cos ot - a sin ot 

k = time index 


T sampling interval 
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tiemouserevation model then is 


ху (К) 
г(К) з [ 1 0 | + v(k) С 
x» (X) 


This model may be written more compactly in vector 


notation: 


х(К+1) = Ф(К) х(К) 


al HOR) X(k) T V (K) 


ur Фа the state transition matrix and H is the measurement 
vector. Тһе on-off keying nature of the signal may be 
accounted for in an intuitive мау simply Бу multiplying the 
observed value, z(k), by the probability that the signal is 
present. This probability is readily obtained from the 
demodulated output of the filter by using the algorithm 
previously derived for mark and space transition probabilities 


for the baseband signal. 


B. FILTER ALGORITHM 

The general filter algorithm presented in Appendix A, 
using the above signal model, was used to filter the (down- 
converted) IF signal. Demodulation was accomplished digitally 
by squaring and averaging the x, state estimate. The 
subroutine used previously for the calculation of Q was then 


used in exactly the same way as for the baseband model except 
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that the indices k and j were advanced by an amount of 

time equal to the delay due to averaging in the demodulation 
process. Additionally it was found that using zero proba- 
bility during space intervals was too low to allow recovery 
when the signal pulse occurred. The probability 0.5 was 
found to be sufficient and was used whenever zero probability 


would normally have been called for. 


C. IMPLEMENTATION AND RESULTS 

This filter was programmed and tested using test signals 
E xbstbcugh +2 dB SNR in a 2 kHz bandwidth, and -15 dB 
ehrouch ll dB in a 2 kHz bandwidth with 100 Hz bandpass 
filtering prior to sampling. In order to determine the 
effectiveness of modifying z(k) by the probability as des- 
cribed above, the filter was also run with the probability 
set to 1. The signal was sampled at 4000 samples per second, 
the (down-converted) carrier frequency was 1000 Hz, and the 
modulating Signal was a square wave with period equivalent 
to a code speed of 25 wpm. Because of the lengthy processing 
time (1 second of data required about 30 seconds of processing 
time), no large sample error rates as such were obtained. 
However, as can be seen from the output, Figures 24 through 
31, the filter performed well on signals of SNR -1 dB and 
above (2 kHz) and -14 dB and above with 100 Hz pre-filtering. 

Figures 24 through 27 are typical examples of the test 


run made using the straightforward filter with no modification 
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of the observed signal. In these figures channel 1 is the 
input signal, channel 2 is the filtered signal, and channel 
3 is the demodulated output. Figures 24 and 25 are the 
results obtained from the signal in the 2 kHz bandwidth; 
Figures 26 and 27 show the results of 100 Hz analog bandpass 
filtering prior to sampling. 

Figures 28 through 31 are the outputs of the processor 
using the modified observation model. Channel 1 is the 
input signal; channel 2 is the input multiplied by the 
probability obtained from the transition probability estima- 
tion algorithm; channel 3 is the filtered signal, and channel 
4 is the demodulated output. The input signals were the 
same as those shown in Figures 24-27. 

Table XIV presents the results of a bit error rate 
analysis made on a sample size of approximately 100 bits. 
Since the test signal was not morse code, no letter error 
rates were obtained. The projected letter error rates 
shown in the table were determined simply by multiplying 
the bit error rate by 10 since this is approximately the 
proportionality factor between these two error rates, as 
can be seen from Tables III-VI. 

Based on this limited data, the conclusion may be tenta- 
tively drawn that a processing scheme employing 100 Hz 
analog filtering followed by discrete optimum linear filtering 


and detection will yield acceptable decoded error rates on 
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2 KHz Bid 


Unmodified Pre-Detection Filter Output, 2 dB SNR, 


FIGURE 24. 
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2 kHz BW 


-] dB SNR, 


Modified Pre-Detection Filter Output, 


FIGURE 29. 
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TABLE XIV 


EEO TATES LOR FRE DETECTION FILTERING 


No 100 Hz BP Filter: 


SNR Bit Error Rate Projected Letter 
C2 kHz) (3) Error Rate 
(ав) (5) 
0 
0 
0 
= 1 10 
Б 2 20 


With 100 Hz BP Filter: 


SNR Bit Error Rate Projected Letter 
(ОКЕН?) . (%) Error Rate 

(ав) (%) 

-11 0 

-12 0 

-13 0 

-14 1 10 

-15 1 10 


NOTE: Rates are based on small sample 
size of approximately 100 


l2 


та ЗаЩ = 





signals of SNR as low as -14 dB in a 2 kHz bandwidth. 
Smoothing of the state estimates at this processing stage 
would probably yield additional gain although an additional 
penalty in processing time would be incurred. 

Although this scheme meets the design objectives 
specified for the processor, it is probably unrealistic to 
assume that such processing could be implemented in real 
time on any existing machine without parallel processing. 
Coding efficiency could undoubtedly be improved but the 
inherent number of multiplications involved in the filter 
algorithm would eventually limit the processing speed. 

One possible alternative is to down-convert the received 
Signal to the 0-100 Hz band, and with a 100 Hz low-pass 
filter in place, sample at 200 samples/sec. With such 


Filtering, the processing described above would be feasible. 
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Vitti. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

The results of the various processing schemes considered 
in this investigation lead to the following definite 
conclusions. 

1. Kalman filtering of the demodulated morse signal 
provides at least a 50$ reduction in letter error rate over 
unprocessed letter error rates of approximately 60$ or less. 

2. Smoothing of the demodulated Signal provides a 
reduction in letter error rate of approximately 30% over 
filtering alone, resulting in an overall Ее of 65% 
over the unprocessed error rate. 

3. Viterbi decoding of the morse characters, using the 
smoothed state estimate for likelihood computation and using 
a first order Markov model of the code to obtain transition 
probabuintwes, provides a further significant reduction in 
error rate if the input error rate is on the order of 10% 
уе. во: торые error rates larger than 10%, Viterbi 
decoding is of little value. 

4. In almost all cases, the processor performed better 
on sloppy code and typical sequences than on perfect code 
and AR sequences. This behavior is to be expected since the 
statistics of the input signal should match those for which 
the processor was designed in order to achieve optimum 


performance. 





5. By using а 100 Hz bandpass filter, pre-detection 
Kalman filtering allows recovery of the code with a resulting 
error rate of about 10% or less for signals whose SNR in a 
2 kHz bandwidth is -14 dB or higher. 

6. The Pickering 230-D decoder (see Appendix C) allows 
recovery of the code with a resulting error rate of about 
10% or less for signals whose SNR in a 2 kHz bandwidth is 
approximately -6 dB or higher, also by using a 100 Hz 
pass filter. 

By making reasonable extrapolations to various other 
processing arrangements based on the results presented here, 
the projected performance of several alternative processors 
may be obtained with a high degree of confidence. The 
processing arrangements presented in Table XV are listed more 
or less in order of increasing complexity and power, showing 
the letter error rates to be expected from each; the projected 
figures are indicated by an asterisk. The rates presented 
assume that the processors, except for the human operators, 
are preceded by a 100 Hz bandpass filter. 

Since no measurements were made for SNR's above -6 dB, 
except for the Pickering 230-D, the projected rates in the 
-4 dB column were all determined by extrapolation of the 
measured values. The projected values for coherent demodula- 
tion (row 6) are based on the presumption that such demodula- 
tion would provide approximately 3 dB of processing gain over 
non-coherent demodulation; the subsequent values for Kalman 


filtering, smoothing, and Viterbi decoding (rows 7-9) were 





PERFORMANCE 


Processor 


1. Non-coherent 
demodulation 


2. (1) followed by 
Kalman Filter 


3. (1) followed by 
smoothing 


4. (3) followed by ` 
Viterbi Decoder 


5. Pickering 230-D 


6. Coherent 
Demodulation 


7. (6) followed by 
Kalman Filter 


8. (6) followed by 
smoothing 


9. (8) followed by 
Viterbi Decoder 


10. Pre-detection 
Kalman Filter 


11, (10) followed by 
Viterbi Decoder 


12. Typical Operator 


13.Good Operator 


TABLE XV 


COEF ALTERNATIVE PROCESSING SCHEMES 


Э 


250 


>50 


>50 


>50 


05 


>50* 


503 


230% 


2@* 


95 


-14 


>50 


>50 


>50 


>50 


>50 


27905 


>50% 


>50% 


>50% 


10 


4% 


2 


2 


SNR (in 2 kHz) 


=! 


>50 


>50 


>50 


>50 


>50 


250% 


>50% 


>50% 


>50* 


5x 


2* 


e 


10 


EIG 


(dB) 


=10 


>50 


>50 


>50 


>50 


~50 


55 


155 


10 * 


4% 


0% 


0% 


Letter error rates (%) 


35 


15 


15 
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3* 


2* 


0% 


0% 


LOS 
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SProjected 


Os 


5% 


3% 


3% 


2% 


1% 


ОХ 


0% 


5* 


O* 





obtained by applying the respective improvement ratios 
actually obtained for non-coherent demodulation. Іп the 
case of the pre-detection Kalman filter (row 10), a 10% error 
rate was actually obtained for the -15 dB entry, and no errors 
were obtained for the -13 dB entry. However, since the small 
sample size tends to decrease the confidence in these figures, 
these rates were increased on a Subjective basis by observing 
the quality of the output signal. Again the subsequent rates 
for Viterbi decoding (row 11) were obtained by applying the 
previously determined improvement ratio. The figures for 
the human operators are the author's estimate based on the 
error rates obtained previously for the amateur radio 
operators. A degradation factor of 3-6 dB should probably 
be added to the SNR's shown for typical field-operation 
performance. 

In summary it may be concluded that the performance of 
a good operator can be approached by a processing scheme 
employing pre-detection 100 Hz bandpass filtering followed 
by discrete Kalman filtering and a Viterbi decoder. The 
performance of a typical operator can be obtained by either 
the 230-D or by coherent demodulation preceded Бу 100 Hz 
bandpass filtering, with further improvement provided by 


Kalman filtering, smoothing and Viterbi decoding. 


в. RECOMMENDATIONS 
A more appropriate model for the demodulated baseband 


Signal would be one which incorporates the exponential 
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rise-times and fall-times of the signal pulses instead of 
the abrupt rise-times and fall-times as used here. A filter 
incorporating such a model would probably have the effect 

of reducing the errors caused by extremely short dot 
insertions during letter-spaces and word-spaces. Better 
estimators for received demodulated signal noise power and 
for character duration mean values and/or probability 
densities should be tested, although it is felt that little 
advantage would be gained from more sophisticated techniques 
in an actual operational environment. A third (or higher) 
order Markov model of the code would undoubtedly show a 
Significant improvement over the first-order model in the 
effectiveness of the Viterbi decoder at sufficiently low 
input error rates. Such Markov modeling, along with improved 
likelihood calculations, deserves further investigation. 

The pre-detection Kalman filter shows the greatest 
advantage, and more investigation of such filtering, and 
possibly smoothing, with the goal of making the processor 
real-time, should be undertaken.  Down-conversion of the IF 
signal to 100 Hz or so followed by a 100 Hz bandpass filter, 


and sampling at 200-400 samples/sec, may be the solution. 
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APPENDIX A 


THE DISCRETE KALMAN FILTER 


The discrete Kalman filter is an extension of the more 
familiar Wiener, minimum mean-Square-error, matched filter 
to nonconstant coefficient multivariable systems with 
nonstationary noise, implemented in sequential, or recursive, 
form. The result of every processing cycle is the current 
estimate of the state under consideration. As each new 
observation is made, the current estimate is updated to 
reflect the information content of this new measurement. 

A brief outline of the derivation of the filter algorithm 
15 presented below. 
The message model is described by the linear vector 


difference equation: 
X(kt+1) = o(k) x(k) + P(k)w(k) 


ее оне поте, or random forcing function, w is a 


zero-mean white-noise process, with covariance 


cov[w(k) ,w(3)] = Q(k) 6, (k-3) 


The observation, or measurement, model is given by the linear 


algebraic equation 
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z(k) = H(k)x(k) + vík) 
where v(k) is also a zero-mean white-noise process, with 
cov[v(k) , v(3)] = R(k)0, (k-3) . 


ПӘРИЕ р СТЕУ 1t 1S assumed that w and масе uncorrelated; 
ено ес упала the initial value of x are uncorrelated. 

The first step in obtaining the estimate at time k is 
to predict ahead from the estimate obtained at time k-1. 
This prediction may be expressed as the conditional 


expectation: 
x(k|k-1) = Elx(k) |z(k-1)] 


Pts nec En. measurement z(k-1) is embedded in the previous 


estimate, x(k-1|k-1), this een may be expressed as 
x(k|k-1) = $x(k-1[k-1) + TE[w(k-1)] . 
Since E{w(k)] = 0, this reduces to 
X(k|k-1) = $x(k-1|k-1) 
The estimate at time k may then be determined by 


considering it to be a summation of this prediction plus 


a correction term employing the measured value: 
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х(К|к) = éx(k-1|k-1) * G(X) [z(k) - z(k|k-1) ] (A-1) 
where 2(к|к-1) is determined as follows: 


z(k|k-1) = E[z(k) |z(k-1)] 


E[(Hx(k) + v(k)) |z(k-1)] 


E[Hx(k) |z(k-1)] + E[v(k) |z(k-1)] 


Under the given assumptions, this reduces to: 


Z(k|k-1) 


|| 
(с 


E[x(k) |z(k-1) ] 


|| 
[x 


$ x(k-1|k-1) 
The filter equation (A-1) then becomes: 
X(k|k) = 9x(k-1|k-1) + G(k) [z(k) - Hóx(k-1|k-1)] 
o Ор и nr l tering is yet to be determined. 
Ее со determine G(k) such that the resulting 
estimates are optimum in a minimum mean-square-error sense. 


Since the error in the estimate is given by: 


КЕ (ЕЮ 


О 





Fm Ar CO matrix OL estimation error is 


V(k|k) ё BL (x(k) = Я (КК) ) (х(к) - #0 |к)) 1) 








Minimizing a scalar quadratic form based on this matrix 
КЕШ Те О ENE opti múm gain с(К). Derivations are available 
in several texts [4], [11], and will not be presented here. 
The optimum gain results in the following algorithm for the 


discrete Kalman filter. 
Gain: 


T 


G(k) = V(k|k-1) H° [H у(к|к-Ин жа" 


Estimation: 


X(k|k) 9» £(k|k-1) + G(k) [z(k) -Hx(k|k-1)] 


Estimation variance: 


V(k|k) = [I - G(k)H] V(k|x-1) 


Prediction variance: 


V(kt1|k) = é V(k|k) ó] + Q(K) 


Prediction: 


x (k+1|k) 


| 
<> 
x 
E 
E 
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APPENDIX B 


BEER SVTTERBT ALGORITHM 


An instructive and informative description of the 
Viterbi algorithm is presented in [8]; the following outline 
is based on this article. The MAP sequence estimation 
problem previously stated in Section VI.A. is formally the 
same as the problem of finding the shortest path through a 
certain graph, called a trellis diagram. In this diagram, 
mus trated in Figure 31, each node corresponds to a distinct 
state of the Markov signal process and each branch represents 
a transition to some new state at the next instant of time. 
In this representation, every possible state sequence 
corresponds to a unique path through the trellis. 


Each branch is assigned the length 


а аа Та (а |р) 


where 


Сұ represents a transition from Ху, ше Ху. 


The P(x,,,|x,) are the sequence transition probabilities 


while Р (2, |5) аге the likelihoods of a particular state. 
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The algorithm for determining the shortest path through 
the trellis diagram iS a version of forward dynamic 


programming and is stated as follows: 
M = number of states; K = length of input sequence 


1. STORAGE: 


k (time index) 


с 1 


| ^ 
|^ 
< 


Io), м (survivor length) 


ООО INITIALIZATION: 


k = 0 


x x )=x x (m cust t mM xX 
x( E о“ x (m) У, Ti 


асада m) = e, m Z x. 


oe) RECURSION: 


алт ссшрщее: 


Bee) Веле), Sor Қ 
b. Find: 
ГО 4) = min ТО А), for each Х 41: 
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(survivor terminating in ху) 





сло оте I(x ) апа the corresponding survivor sequence 
~ 


ри 
Е 


Ч ЕС К EO ktl and repeat until k = K. 


As an example [8], the recursive determination of the 
shortest path through the trellis shown in Figure 32 is 


shown in Figure 33. 





k=0 k=] К=2 К=3 К=4 К=5 





Figure 32. Trellis Diagram with Assigned Lengths 
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4 1 k=] 
1 


ў 


3 
k=4 
3 
4 
К=5 


Pace tos. Example of V1iterbi Algorithm 


F à 








k=3 








APPENDIX C 


PICKERING 230-D PERFORMANCE EVALUATION 


The Pickering 230-D is a currently available automatic 
decoder and transcriber which includes some front-end 
processing. The front-end consists of what is essentially 
a type of coherent detection scheme with a bandwidth of 
about 180 Hz; the center frequency of the unit tested was 
ас 875 Hz. The letter error rate versus SNR of the 230-D 
was determined in the laboratory using the model KB-1 
keyboard to send perfect code. The error rates were 
determined for both a 2 kHz bandwidth and a 100 Hz bandwidth 
and are shown in Table XVI. A performance evaluation of the 
230-D for error rates obtained using actual operator- 


ЕРИСИ Е Ееа code with no noise is available in Ref. [2]. 
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TABLE XVI 


Pickering 230-D Error Rates 


SNR 35 wpm 30 wpm 25 wpm 
(in 2 KHz) Е 100 Н2 2 kHz 100 Hz 2 kHz 100 Hz 
(ав) (3) (9) ($) (9) ($) (5) 
6 1 3 1 0 0 0 
0.4 4 ES 4 6 0 0 
0.9 29 Be 22 13 21 6 
-2.0 = Г ES 11 - 6 
-4.4 - 7 = 7 - 4 
-6.4 = 15 = 15 - 5 
55-9 = 21 = 27 = 22 
E05 = 60 - 78 = 61 
-12.4 - _ - - - - 
Notes: 


1) Sample size: Approximately 100 letters in each case 


2) Perfect Code 
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1; 


هم حو هم 
г н» о 00 У W‏ 


حم هم 
ғ (9‏ 


B5: 
16: 
7: 
ше: 
шо; 
ро: 
21: 
BO 


>» N rO 
1 4 а) 


NUN & с Е ЕЕ с Е Е Е ОЕ 600 (а) ОЎ (а) CO CO Y Y Y ба) ГО ГО ГО ГО 
— Су 540 20 ~ СЛ 4 ӘӘ ШО О +J O" U) £ Q) (N — O Оо УСО 


ee зе ее ее ее 40 Ae ве ee ее зе ее ее ее зе зе ее ее за зе Фе 49 Фе ее os oe ее 


г 





COMPUTER PROGRAMS 


AORTRAN [5,69 


ща се сесессссссссссссссссссесесессесссевсссбвровеоаас 
С С 
С THIS PRÉGRAM IMPLEMENTS THE POBST«^ETCCTISN KAMAN C 
C DUERCOAND LDIMCAR SWRGTHE3, ALENG AITH THE CEMED € 
Е SIGNAL PARAMETER ESTI“AT92 ALGORITHMS. THE FLIM С 
С UTILITY SURBRSUTINES ARE USED: H 
c ӘЦЕРЕЗІМ = READS DATA ON INPUT TAPE. C 
С DAL = O/A CANVERSISN AND SUTPUT REUTINE С 
С THE F8LL8WING SUSSBJTINES ARE PRAVIDED: С 
C BACK = BACKWARD FILTER REUTINE С 
С STATS= ESTIMATES SIGNAL PARAMETERS С 
С VARW = COMPUTES THE VARIANCE OF W(K) Ф 
С С 
а ШИ БАГАНЕТЕКС АТ ТІРЕ ВЕ EXECUTION ARE :! С 
G INITIAL ESTIMATE S&F Т С 
E THE PARAMETERS T1 AND T2 С 
С ВЕР ТО ВАТТ DENSITY PARAMETER, (Dal C 
C THE DASH*»DURATI8N DENSITY PARAMETER, DOAS4 C 
C C 
NNNNMCICCCCECEECCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCESOSOE 
C 


РОА ЛАО РИТУ МЕ АМУР МЕКА МЕ А XD, VNXTEMNXOMERTT 
ВИМЕ МО РОМ XYHSMNI (50), XHSM(159) Х45 УУ 2(120),Хх45м2 (102) 
По EVARP (750) a XMATE € (250) se Vener (ЕВО) 
Ра ХВ! 25)» РУАРЗе (125), ХЦАТ 21125) 

DIM ES ЕЕ (ООО), ХНАТР (250) 

ВОО АТВ (ОБО )„РМАКВ И (250) Уа (е5 Сте 
DIMENSTON МЕБІТСЕ; 2,123, ТХАМӘ(11) 

ПО КО а ЕК /XNYT yXNTC XHSM KHSTBZ T JH 1157 J27 7 тоо 
DISSI 972 7 DNS; NSHÖORT 

(ПАН АРСИ МЕ ЕМ, УЕВ 11 

(ОПАО РАВНО СЕЗАТХАМБ 


МАМЕ 157 LOSTsDDASH,TsTIsTe 
E Sree teansiliGN PROBABILITIES FER VITERB] OECOJEN 


TRANS (4) =*ALEG(0+682) 
TRANS (2) =-AL9G (00318) 
TRANG (3) "AL 8G (0504) 
TRANS (4) ==ALNG(0+396) 
TRANS (5) =-ALOG(0025) 
ТЗАМС (6) = - А65 (0605) 
TRANS(7) s -AL8G(04 25) 
BE сос 251 
ТВАМС (9) <-А| 26 (С „ба ) 
TRANS (10) =-ALOG(0+316) 
TRANS (11) ="AL9G (O01) 

39 CONTINUE 


Є 
8 INPUT PARAMETERS 
С 
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92: 
53: 
58. 
DE 
56б: 
EE 
Son 
БОР 
6б; 
61: 
Бе: 
63: 
6%: 
6S: 
66: 
G7: 
62: 
Go 
70: 
21: 
721 
л 
74: 
zo 
76: 
TA e 
pu 
79: 
80: 


82; 
oon 
84% 
on. 
86! 
87! 
58: 
49; 
30; 
S 
92% 
93; 
24: 
J5: 
EOS 
27: 
7, 
л: 
100; 
TOIS 
102: 
198; 
104: 


ШО ТТЕМТЕК Т) ТІ, АЧ02 Te (MSEC)! 
а ФЕЦТЕК LOGT AND DDASM (MSEC)' 
ИИ 


ШЕ НЕ VARIABLES AND INDICES 


СУ СУ СУ 


TsT/2* ;ј LDOTsLDB8T/2« 3 DDASH=DDASH/2: 
IM uc" Т2-Т2/2. 
A*sLDOT/T 
ВерГРОАЗН/ (Зе т) 
аа МЕАЧХВЕО.Б 3 MEANX=0¢5 3 XT=T ; 11=7 
= т 
IX=0 
KS1=KS2=0 
ИЕ МААЕ] e 
51 +6р: .5 
5512552=44 
МАХ ] =ММХ2=0.5 
IM1=-]M2=151>-152=0 
АЕТ РЕ ротора 10 
У = а ІО 
NELEMENSHART =O 
ЈАЧИ <= ево = | | 12] 1220 
SMN]=S5MN2=0+ 
XHSM1=XHSME=C « 
ЖОС Л Л О ; XHSMNe2(1)=08 
[хМ=0 
PVAR=EVAR=VARY=0D005 
КОШО ШО ЛОК УОК 150950) JINDF=0 
Gates (ИДІР-Об ; |хн| #0 
muse ; XT1=XT2Ə=T 
XTD1=XTD2=3*x*T 
XILISXICAST 
ХТА=ХТСЕТ 
DA 2 К1:1,% 
_ Peewee І], 2 

OE 2 КЗЕМ,12 

2 ЕЕС ККЗ у= 0 

38 CONTINUE 


C 
C INPUT DATA FROM ТАРЕ 
C 


112 CALL BUFFERI*':Cí,1,IBUF^, 500, IERR) 
ша аа аба ў) 49 Та 1 
ЕТЕУ 1,210,112), TERR 
IIIE SINCE Hise 
О 1 =( КОС 1 зе 501 
INDE2= INDE*250 


E EGIT TESCESSING 


ВВ ТРЕХ = 1 МЕ! , МОЕ? 
ООВ Е С] ЧОЕХ ) )/ёж+23 
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O5: 
106: 
TO 
108: 
109% 
TO 
LI 
112: 
Ig: 
S 
015: 
es: 
фиг" 
118: 
Ic 
129, 
el: 
с: 
В" 
124 
125: 
126: 
127: 
122.: 
feo. 
180: 
зү: 
132: 
SE 
134: 
195, 
136: 
137%, 
135: 
Mc 
140: 
141: 
Duo! 
143: 
144: 
145: 
146; 
paz. 
148: 
149: 
Den 
Б 
158: 
Це зе 
153: 
Do 
156: 
157 


OOO (3693900 


Cy OVC OG 


OONA 


E 
С 
С 


Dein 


СҮҮС) СУГУ 


Шар ОЕ IUTRUT; X9UTI 
ХӘОТ1:0. 

12 (Х1. СЕ. МЕАМХ) ХВОТ1=1е 
NS ERAS A ME INDEX 


JAS 1X+1 
cines 0). 1X=1 

ISS MES TATNG ALGGRITHY UNTIL SUFFICIENT 
DS NEP (500 SAMPLES? 


| ХКЕД ХМ + 1 

IF CIXNeGT e500) ЈЕ РО 

STBRE INPUT SAMPLES IN REVERSE GRDER FER BACKWARD 
FILTERING 


¢ 


ІЕ(ІХ»67.125) Х81-ХВВ(1Х-125) ; 
XB12XR(1X) 

ХВВ (126-1Х) =ХР (251-1Х) 

ХВ (251-1Х)=Х1 

CALL STATS 

VARXF=VARX, 


GO TO 9 


FORWARD FILTER ALGORITM 
GAINFEPVAR/(PVAR+VARXF ) 
ЕУАЗг(1»-САТЧҒ) ҰРУАЗ 
XHAT=XHATP+GAINF *(X1-MEARX=XHATP) 


THRESHOLDED FILTERED SUTPUT, xGuT2 


ХЗУТ2=0. 

IF CXHAT«GE • ОФ) хай аз І е 

[ХНАТ= = 1 

|Е ( ХНАТ „ЗЕ O, ) IXHAT=1 

KINDEKIVOF 3 JINDEJINDF ¿ IXHL=IXHLE 
Т=ХТ 5 TD=XTD ¢ LOBTSA*XT 3 DDASH=B*XTD 
ХНЕХНАТ 

CALL VARw 


ОРз(к (РУБАКМЕКР) 
РУАВ + ЕМАВ + СЕ 
KINDF=KIND 3 
ХНАТР=ХР 


JINDFEYIND ; IXHUF£IXHL 


FOR 


STORE FORWARD-FILTERED ESTIMATES ANO VARIANCES 
В Б 


ЕМАЗЕ | ЕЕУАКРФ (| Х); 
EVARFOTX) SEVAR 


СУАНГО( ІХ) = ЕМАКЕ (ІХ) 
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Су СУ СУ СУ 


ЭЕ 


493 


OOO 


EX c0 


XHATF1=XHATF2(1X)3 


ХНАТР(]Х)ЕХНАТ 
CALL BACK 

CALL VARA 
QBtQe(FMEANSe*2) 
RINCIFKIND > 


ХНАТР 2 ( | Х) + ХНАТЕ ( | Х) 


JINDB=JIND 


e 
? 


IXHLB=IXHL 


REVERSE STERED VALUES FRA" BACKWARD FILTER FOR 


SMBQTHING 


IF(1X«GT.125) 
ЖАО Хо 


ХНАТВЕ (126+ | Х)} =ХНАТЕ1 (251-1Х) 


CoU GENUE 


XHATB1(251-1X)=XHATEP 


ГЕ ХЕ 514122) 
РУН | «РУАЧВ1(1Х) 


РУАЦВЕ (1265 | Х) <РУАКВ!| (251 - IX) 


CONTINUE 


РУАКА 1 (25 1. | Х) ЕРУАКВ 


ХНАТВР + ХР 
PVARB=EEVARB+OIB 


СИЗИ ING ALGIRITAM 


AKSEVARF1/(1.+EVARF1/PVB1) 
EVARS=((1o-dK/PVB1)**2)*EVARF1 + 
ABK=(XHTB1/PVB1)*EVARS 


ЕШРІ-ХНАТВОСТХЯ 1295) 


РУВ1ЕРУАКЗЕ( | Х- 125) 


ТЕ (ТХ. ЕСегь0) WEK=s0-6 
ХНАТЗЕХНАТЕ 1 / (1. +ЕМАКЕ | /РУВЪ) + В« 
ОО) XHATS=S02 3 TS=XT ; 


AVERAGE AND 


BEER 7555716 
ЕТЕ! 
IF(111+GTe150) 
Вт) 

I | рат. 50). 115 
Beer») IT) ГІ 
м1) 
МБ (ІТ) =ХРАТ5 


1 


1 


111=1 


Ato 1=XHSMI+(XHATSSXHS1)/(0+75*TS) 


У Е ХНОМЪ 
= о. 67 705 
112=112+1 
тест» тоа) 
НЕ Сте 1 Те) 
INE Тат 
Ш е» 100) 
XHS2=XHSMNel lle) 
ин | $ XHATS 


1122] 
I 


113=1 


1,389 


но" Тоша бе 


аа 28 


(иКнжё) /РУВ] 


ОБЕХТО 


s ЛОТ ЕСТІУАТЕ AT EACH POINT PER 
IS PI T R ІМӨсСІКЕІЗЕОеО СОӘМРОТАТ!19М65 





с: 
pO 
Boe 
214; 
дей 
l6; 
= 17: 
ale: 
219: 
ee): 
Bon: 
222! 
Bes 
224: 
газ" 
226: 
227. 
гга; 
225: 
230: 


Sl 


232, 
сс: 
234: 
299. 
236: 
cuo rt 
238, 
сло; 
040: 
241; 
puc. 
сүз: 
244; 
Bus: 
246: 
247; 
248: 
249: 
250: 
gr: 
268: 
со за 
азі: 
go 
CS6: 
ET 
сол, 
2509" 
COD 
eed 
262: 
Eons 


ООО С С СУ 


J20 


сащ 


XHSM2= XHSM2+(XHATS-XHS2)/(0.67*T1DS) 
AS PDT ISO AS ME 

ЈЕНИ ТО 

ПЕ ЕЕЕ ГУСЬ 

IJisII1*5-IIT 
ТЕСТ ЕО) 
поје - Тео 

ПЕ O) 
XSuM=0 

LO 5] so; 10 
XSUMSXSUNM*XHTGCI) 

к (І) 
XSUMSXSUM+XHATS 
ЖЕЗ нА ПО 

X3UT=08 
IF(XSUM*GT+.D0:*0000001) 


Е 


шанае 


x9uUT:le 
ESTIMATE T AND TD 


XNT IS THE MARK DURATION COUNT 
Reel ET SPALE DURATION COUNT 


KND=0 

ПС РА и ТЕ ) • „Се - 0,39) 
ТЕ жест СЕ е 0.99) 
МТЕО 

ӨБ 7329290 

ХМТЕМТ 

NT=0 

IF COXNT«LTeO4S*T1)ePRe XNTeGTe6**T2)) 
о ЕТ, + те/е«) G9 TB 991 
JTDsJTO-*1 

IA o.O)JID2s0 ; 
ПРО Тре јео) JTD=1 
XTD1=XTD1+(XNT-XT201)/JT10 
аца а? 

JC JO O, [QL] 10) JITOÍ=0 ; 
ПЕК трена ўе 1то2=1 
Е хрретохмтехтоеуи тог 
XTD=XTD1 

пары ас те р іа ХТОЕХТОС 

IF CXTDeL Te (2e*XT)) xXTD=2+*XT; 
а Је трев=о 


KNDO=1 ; 58 
ЕИ со те 990 


ТО 


IO а 


хХтрагхтој 


xps е 


xTD1=:XTD2=XTD: 


Comes o 

Sut 

ТЕА БЕО Јтезго ; XT2sXT1 
Неле се 20) Jys] 

кача LAN TX TIT 
JT2=JT2+1 

Имате ее Ој отео ; XT1sXTO 
JI Мое .20) тег) 


XT8223XT2* (XNT»XT2)/JT? 
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06%; 
265: 
сев" 
267: 
БЫЗ: 
269: 
DA 
ЕЙ: 
ЛЕ; 
t 
DIEN 
aps. 
4 
r 77: 
et 
279: 
END 
Pd t 
202 
сая: 
08%; 
es^ 
CS6? 
ES e 
са! 
259: 
zen: 
po 
22 
газ. 
294: 
p 
= 26: 
соу: 
298: 
EE 
3005 
201: 
202; 
303: 
S04: 
325: 
306: 
o e 
SOR: 
309: 
SOn 
ӘМ: 
зае: 
виза 
31%: 
1\5: 
eve. 


390 


100 


ШЕ 


(C) CY C) 


ee 
J 


Do 


GON 


EIS 


SI 


OOM 


(OOOO 


ХТАЕХТ | 
Јали о Пе ХТАЕХТЕ 
CONTINUE 
ОИУ 6E«OV99)-XNOs1 ; 69 T9 1000 
AA INPICENTCKI 7 350 T8 1090 
ана акад 1090 

O XNTC=NTC 
NTC#9 
ПОРО ТТИ ова СХМТС.е ОТ .ебекто2) ) 69 T9 1099 
J JUD е 2122۰ JIC=JTC+1; G8 T9 1091 
G921792 1020 

Вто) JTCe=O ; XTC2sXTCI 
E ооу JTC=1 
XN NIU XNTCTXTCIZJTC 
Све јеж! 
ПОШ о) 0112-50 ;) XTC1=XTC2 
Пе. 020) 21С2-1 
Oe OC XNDCEXTICeY/JTCO 
ХТС=ХТС1 
ЕТ. ТС} ХТС=ХТС2 


РЕ ОККО ОЕП ТО SELECT VITERBI DECODER 
OOS TINS PROGRAM, EXECUTION 


DENM SE SWITCH 1) 109251093 

à IF(NTCeGTe*e«*TS) KNDs1 i XNTCsNTC ; NTCsO 
ВЕ АЕ О) CALL LIKELIKHIBD 

ВС E 
МЕТА» ТС) ИС. 
XSUTLEXSUT 


ШИЕ АЕТГЕЗЕЕС ТО ОХА REUTINE FSR ANALOG RECSRIDINS 


Nele =e xX 17100. 
NI ТТ ХӨПТО,ХНАТ»ХНІВІ, ХБАТ5» ХӘШТ;УАХХ;Ҙ 
1GF »GAINF, TNSRM) 
CONTINUE 
60 18 38 
ОШОО (ООУ Емо SF RUN ; HIT = TO GO ” 
ЈАЈА 0101) 
GE те 39 


SUBROUTINE BACK 
ШЕ КАО Т МЕ IS THE BACKWARD FILTER 


INITIALIZE INDICES FSR BACKWARD FILTER BASED SN 
BEE SIDE FERRWÄARD=-ETILTERED ESTIYATES 


ЕТ. Е.) 69 TO 490 
ЦК ТАРЕ ЕС)" « [ора=0 7 58 T9 360 
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3171 KINDB=XT=-KINDF 


SOS ЕТО) KIND3=XTO=KINDF 

sy: РАНА T.00) KINDE=O 

320: 360 ТРЕЕ lp) JINDG=0 ; G6 TO 361 

Sal: JINDA=XT+=JINDF 

322: Pee ote) JI NDC =x TO=JINDOF 

sen. IFIJINDB«LTeO) JINDE=O 

pue 361 CONTINUE 

oos Xu Wo NIE IXHLS5zIXHAT 

СЕ; PVARB=EVAR 5 VARXB=VARX ) МГАМХЗ=МЕАКХ ; СВЕСЕ 
ce ТЕ с TOR ; PrEANS=SFMEANB 

T t FMEANBEFMEAN ; TB=XT ; TD3=XTD 

329: 490  CENTINUE 

p C 

Ew C FILTER ALGARI THM 

БЕР С 

332, GAINB=PVARB/LFVAR3+VARX?) 

334: ЕУАНВЕ (1, -ОЛТУВ) “РМАКВ 

So ХНАТЕ + ХНАТВР + СА | МЕ « (XS1-MEANXS=-XHATSP) 
336: JXHATBS+1 

moy ІР (ХҰНАТЗ СЕС» J) [IXHAT3Ə3= 1 

332: | ХНАТ+ | ХНАТЗ 

339: В как + JIND=JINDS + IXHL=IXHLB 

340; ош ел СБ + T=TB / TD=TDRB J LO9ST=SA*T 2 DDASrHsS*T23 
341: ХНЕХНАТВ 

342: 450 RETURN 

Emo С 

3441 C 

ее С 

346: SUBRSUTINE STATS 

247: c 

348: С l CUTIE CALCULATES THE INPUT SIGNAL MEAN? 
349: С VARIANCE, AND ESTIMATES THE SIGNAL AMPLITUDE AND SNR 
BSO: С 

ӘЗІЗ С 

ЕС INPUT MEAN AND VARIANCE CALCULATIONS 

post C 

Etpe = МТ 

8565 IF (IM12E9-590) 1м2=0 ; МАРЕУА! ; МУХРЕМУХ1 
ESO: шае „Епо ВОС) 1411 

357: MNX1=MNX1+(X1-MNXx1)/I]M1 

SE: Via noe Ot XP ME ANX VaxceVAL)/IMI1 

EO ies i eri 

360: 222222500) Г/11-7 ; VALEVAe § MNX1=MNXC 
361: ие. Ое 1000) [Meri 

362: ЛК + СХ|-"имхеј)/1Ме 

363: МАСЕУАЯР + (| (Х1-МЕАЧХ) »#О "МАР ) /И | М2 

36%: MEANXEMNX 1 

365: МАЕУА + 

366: Те 11) МЕАЧХЕЧМХО 3 VA=VAC 

3673C 

355: С SIGNAL AMPLITODE ESTIMATOR 

SE C 
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3703 
37 1 
а/о 
уз 
374% 
3793: 
ӘЛЕ” 
de 
378: 
B72% 
220: 
381; 
вази 
363: 
3843 
gas: 
356: 
E 
388: 
389; 
СОВЕ 
Ec 
зэв: 
sce: 
DA 
DEL 
3236: 
2397, 
Sam: 
295% 
>00: 
401; 
402: 
&03; 
404; 
405: 
806: 
407: 
4OR; 
409: 
410: 
al: 
aoo 
412: 
414: 
415: 
416: 
417: 
418: 
O 
ipeum 
421: 
нее. 


MAND ODO 


CY yy CY CY Cy O 


51 


10 


IF(X1+LE+MEANX) (8 TO 80 
psu 

ЈИ | син SIO! |52-С ; 52-51 
БОО Т КК ООС ISIS 

s Xo) ) 7 1S) 

ВЕ о! 

паа гама дс д 165180 ; 51550 
Laer О0 0) [5251] 
ӘСШФЕЧ(Х1-52)/1652 

ок» | 

и се о Пе 151) 5-52 

| 

CONTINUE 

СЕЕ КЕТТІ 

ПЕС 995500) КО2е0 4 5525551 
ЈЕ (К51•20•1000) KS1s] 

Be soot (41 =SSl)/KSi 
(асака 

шаша аа 05500) «51-2 1 5515557 
nass 50:109000) KS2=1 

55-5 -25(Х1-5502)2/Кк52 

А. 

J (O СТЖК51) 555552 
CONTINUE 

FMEAN=S-SS 

И АМЕА ОТРУТУ) ЕМЕАМ=ЕМІМ 


Ne ISE VARIANCE ESTIMATOR 
УАКХ=\МА-0 › Эх (ЕМЕА\Мжжёе ) Иф. 


ТР (МАКХ | ЕО.) УАКХ=0. 0000001 
RETURN 


SUBRBUTINE VARW 


Q ESTIMATION ALGORITHM 


р = ОХ ка FOR MARK PROBABILITIES 
Vu TIAE NOE Xs? Pas SPACE PKOSABILITIES 


ПЕ aO FILTER? INITIALIZE KINOsJINO Та 60) 


ХРЕХН; РЕ: 1.0 
ЈАЧИ ЕО) GA 70 10 


ЈОАНА ар.) СТО Та АМОС 1 МО От • те 28711) 


ЕО Та ще 
РОС | ХНАТ- 1 ХЧ|) „СТ СО) ТХНАТЕ- | ХЧАТ 


раса пе Хе ето ОС) е АмО у(кТ МО. СТет- 1 ОТ) ) 


1 | ред; GO TO 10 
ІР С(ІХНАТ-І ХЫ) СТО) | ХНАТе | ХНАТ 
рам eAT eile 1) KINDEKIND4+1 
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° 
г 


° 
U 


XP == XH 


xP=-XH 


01 1020; 


5956: 





ips 
42%: 
425: 
426: 
нед. 
408; 
го; 
+30: 
з 
432: 
433% 
434: 
435; 
4 361 
437; 
432: 
439: 
440; 
451: 
442; 
443: 
444: 
час; 
446; 
447: 
443: 
149: 
150: 
iol: 
ict 
453: 
454: 
155: 
456: 
457; 
Нов: 
52: 
4603 
561: 
462: 
463: 
464: 
„65% 
466; 
467; 
4683 
469; 
470: 
ut 
inl ct 
4731 
474: 
4753 


С 
С 
С 


ms 


ще 


Ша 


ша 


17 


18 


20 


Е 


Е > 


100 


РАР МОЈИ 3 GO T9 20 


MARRAS ELITES 


ХКЕКТКО 

Mae LOST) 
Gen io ic 

G=0 

IF CKINDeEQe1) 
(АВ Т8 100 
NEXU E +L OUT) ) 2 ANDe (XKeGTe(T=LDOT))) GS TO 13 
GO T8 1^ 

РЯКО 1 - (Те ХК)/ (цех! ОЗТ) + 0 75 

РЫҚЫ15:1» - РИКО01 

Ree ОТР РЕРЈКОЈ, G9 T9 100 

IF-CUXK LE «CTDeEDDASH))«ANDe CXK*GE« CT4LDBTD 2) 
НТО 100 

Ее СТО родн Ам 
G9 T9 18 
PUKO1=(TD-XK)/(2.*DDASH) + 
РАКК121» "РЧХС1 

(= РукОј хХРАКМ]; РзРИКС1: 
ЕЕ: 29:76. 100 


(G да 


ХРЕРМРАЧ/О, 


250.5; 


«(XX «GE«(TDeDDASH))) 


59 
O 5 


Ge. ТО) 100 


САРУ ВАВ ТТ1ТЕ5 


XJ*JINO ; РзОеб 

ee wel wt |e ВТ У) 
по њег 

020. 
ПРУ из оте Бе Сре) 
(НОО 

ДА Уто те рату) «АМО (ХЈ•ебЕ • (т=- рет))) 
50 78 24 

РякоОоОт (+ Х.)/ (Рех| ОТ) + Ое5к( | е + Т/| 26Т)) » | Ре / | етъе/ 17. 
РикУ0Ов |. -РУКОО 

GEPWKOO*ePWKNOS P=ePWwKO0N+0e53 08 ТЗ 100 

IF (CxXxJeLEe CTO =DDASH) ) eANDe (XU*GEe (T+LO9T))) 
Pow то 100 

ОТЕ ПАЗ) «АЧ (XJ GE. (TD-DDASH?))) G$ 
GO T3 28 

РикКООз ( («ХИ (2. «БоОАЗН } +0,5* (1 .+ТОЭ/ЭПРАЗН) ) } *О. 4+0» 2 
РУКМОг ] •• РМКОС 

Ог-РАКООяРАҚЫО Ргрчік00%0.5) G8 
Ia) AOS ; “639 Te 130 
асада аа тэа veo 
РУКМО= 1. «РикКОС; QEPnKOO#PAKNOS 


502 el 


XP="FMEAN/2 e 


СО ШЕ? Es 


irs. 


а 


№ 100 


P=PWKOO+005 


CSNTINUE 

МИ В О ВКО: %-0.25 
От) 115220 ) 220.25 
ІР(РЬСТ, |.) Рк1.0 


]ХНЕ = | ХНАТ 
RETURN 
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каа 





MD ODE COMPILATION 


: аша сата е Сас НЭ р 
= (С 
m В ИЕ СА СЛ АТЕЬ THE LIKELTHSPDS ӨР EACH 
“с CHARACTER „ USING THE METHOD GIVEN IN THE TEXT GF 
№. С THE THESIS 
6: E 
pt «С ESI IS THE DOT LIKELIHESO 
ва Е ФГ) Іс ТЕР EVEMENT=S°ACE LIKELIHIED 
E С ЖИІ ІС ТӨР ІСТТЕК-ӘРЛСЕ LIKELIH95D FIR 
ше С СО РОГАТ ТВ НО MADE ON A CHARACTER VS 
и. С CHARACTER BASIS 
в С маса сс шас "ГЕ BE TTEX©SPACE LIKELTIHOSD FUR 
үз; С COMPUTATI@NS INVOLVING THREE SHORT 
ша, С ‘ CHARACTERS 
ШЕ 8 МЕКІТ(5,1;,М) IS THE CHAR VS CHAR DASH LIKELIHESD 
16: C MERIT(4s20N) IS THE THREE SHORT CHAR DASH LIKELIHSS2 
ше С 
ща С 
ща БЕЛІГІ МЕСТТ 
20: BEE IT OSEMERTTLSSZIL2DIIXFERT (3) 
ey IDE ЗОО МКП БО), ХНОМОрЕ 100) АЕБТОРЕ (8) 
с; ОИСИ NOMASKCLTY)SNSPACEL11) 
ee: асада TIZXBNIGQXNTCXHSM,xHSMD,,IJ»IJi,»,IJe,IJEb541J«, 
С»: I SE ANS INSHT 
2t адама ТСКОЛЬЕГЕУ;МЕРІТ 
26: COMMON/BLECK4/MARK se NSPACEs NCHAR 
mt NCHARZNCHAR-*1 
ec: MARK((NCHAR) 2e *XNT 
29: NSPACE (NCHAR)=2+*XNTC 
ӨШ? С 
EIS ( MS ТЕМТАТІУМЕ DECISIONS BASED ӨМ THRESHSLDED 
Bes С ама HER OUTPUT 
BHO С 
34: Г. T.O 55) ОВ ТО 103 
35: IF OXNT eGEelee TS) FOuR=XNT;j NSHORT=0; 2Е5ЈМ20•; 
36; ШЕСТО 102 
ay: Пеи і Тео т) АКБ, (ХМТ СТ. О.Э) ) EDUR=XNT; 
38; ПЕСЕН 101 
а: 103 IF(XNTC*GEec«*TS) EDUREZXNTC; NSH8RT*0; ESUM=O08s; 
40; в о 
41; IO eet) ee eOR1S) «ANDe (XNIGaGle0+5)) EDUR=EXNTC 
42: C 
#3: C COMPUTATIONS FAR SHORT CHARACTERS 
45: С 
45: 101 NELEMSNELEM+1 
46; NSHERTENSHSRT +4 
47: ЈЕ (монант«ецат«еа) NSHAIRT=0 
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ч 
U) 
CY C» C C3 


Sa 
86: 
87; 
ва; 
89: 
90: 
SIS 
ӨЛЕ; 
93: 
ОС 
SEC 
г 
Ы 102 
аа 


ОБН 
100% 


ESTARE (NSHORT )=EDUR 

E SUMBIE SUM +EDUR 

ISA SO ESUI= ESUM=ESTO RE (NSHERT-3) 
TM1i=10e 

ПТ вое > CEDUR+TS)/(2.*TS) 

Тео. 

N ра 9 =1 «= (725-2008 )/(0.67*Т55) 
ХЕ ХЕБ Л ИЕ МЕЖНЬ+Ов 5 } *ТМ1 * (1.-ТР2) 

IF CXM4ERIT«LT»O.) ХМЕК!Т=0. 

о Тов) XMENRITs1* 

МЕ АНА АЕ ЕУ) ХУЕНТТ 

МА ТАИ НО ОУ Јо ЛЕНЕ до 0 + 5 ) ) ТИ] + (]• »Тма) 
ТЕУ КЕЕ Т» 1.0») ХУЕКТТ 50. 

ІРХЕКІТ. СОСТ. 1.) XMERIT= 1o 

МЕКТТ (2, 1, НЕ | ЕМ) + ХУСАТТ 

AS IS) МАК ТО (1 ее ТМ 14 (14 -ТМО) ) 
Td) XMERTT=0O00 

[0264 EG) ХМЕКІТЗ е 

MERA NELE A) S MBR IT | 
A ава А О О ЈО ГРУБА + 065) 1 ТМ] » (Те "тма 1) 
ТЕГ. О.) ХМЕН]1Т=0. 

ТЕ Еро пе С Тате) ХМЕРТ ТЕ | е 

Kern NELEN)=SXMERIT 


МАВАРА ШОККО seo) COMPUTATIONS FOR 3 OR MORE 
ЗОО А се Б IN SUCCESSION 


От Те =) ша те зоо 

и РЕМ Те 3) се TG 300 

DOES Dm] Э 

КЕ О О МАХОМ ТО Ју ДЕЦЕМ- 3+ 1), 


3+1 )) 


XMRTs(1«-XMERT(1))*(1« "-X^ERT(2O))*(1«-XV"ERT(3)) 
XMRTEXMATxAD0 33333 

ІМ2г1. 

рашае аец Тет 05) TMeziwee-tTOS-ESUM)/(O..67*TDS) 
DD Теје) ТЧ222. 
XMERITz(XHSMD(UIJ2)/FMEANS*O» 5) * TM2*XMRT 

JF CXMERITeLTeO«*) XMERITsOC« 

IF (XMERT T«GTeie) XMERITs1* 

ПӘКШЕСТЗе, ХЕГЕМес)-УМЕЕТТ 
XMERITzs(1«-UOXHSMDUI J8 ) /FFEANS*0«5) ) * THO * XMRT 
DESC Иве] те ТО) RMS O. 

ПР ео Теје) ХМЕБ172=1• 
ЕРИ ОУ МЕГ Емо) = ХМЕВІТ 

68 Т8 300 


ТА ТН ЕВЕ LONG CHARACTERS 


ШЕШЕМ ЕКЕ ЕМ+З 
IJK21J9-EDUR/3*« 

E VEO) TJs1590«1JK 
Јоде СО)” З е 
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arm cy 


od 


la). allen eng Nee 
те 


IA сына 1s) пир њис • = СЕрОК+-Т5О) /(2• #15) 


ТИСТІ), ; ІГ(ТУ1. 57.1») 


ТМе=1. 


те 


ЈА А о вео пе = = Страст) ССеб7%105) 
ИРЕ ЛС "Е Ао" Се о у ктманк Тае ТМ) 
|Е (ХУЕР Те Тебе) XMERITEOS 
ТЕРЕГІ 14) ХУЕКІТЗ1» 


E ООМО 2) -Х“ЕКІТ 


X NEN cCMPOLISSEYX/FMEANS*Oe«5))«TMOv-(l1e- T^i) 


Г.С.) XMERITsOs 
ТЕ SI sO. s) XMERIT=10+ 


Е селе } = ХМЕКТТ 


ШЕП жена ГАЗ ЭК FOR DST AND ELEMENT SPACE 


БЕ TODS 


J G ми Зее! Тебеъ“Т5) 9" TO 3CO 


ХМЕРТТ=(ХН$М (1 }ИРМЕАМ9 +05) * (1. ТМ (16. м1) ) 


КОО Гре Т.б.) ХЧЕКІТ:0: 
ШОСТ о у ХЧЕБІ1г1» 


ШЕГУ, ЕС сео југ ХМЕНЕТ 


КОДА ж ХНСУСТ | )ИЕЧЕАМЪТ СО) ) + (| е" ТМОХ (Те ТУ 1) ) 


Е ( ХМЕНТТ.! Т. Се) XMERITsO« 
Шаша Геба Те е) XMERI Tele. 


РАВАН о МЕЦЕ Мер уз УМЕК]Т 


ЕВЛИ ИСО Је ИЕКРАМО + Сер) ж (1 ее Ток (1е ТУ 1) ) 


ШОО КК Тор Тебе) XMERIT*O?* 
ІА КАР ІТ «ОТ ае) XMERIT=10+ 


МӘ ре ро Ту НЕРЕМе Г УЕХМЕК ТТ 


ХМЕР | Тя (1ее (ХНОМ (Тук) ЕМЕА МО + 0 +Б)) # (1 ТМ2ж (1е-ТУ1)) 


РТ Г. Т.О.) ХМЕБВ]Т=0. 
ПАДАЛА Пе а пе Шеј XMERIT=1e 


(а по ме См у = ХМЕНТТ 


ХМЕВ1Т= (ХН$М (1.1) ИРА +0. 5) * (1. -ТМРх (1. =ТМ1 ) ) 


BREXONEERITeLT«Os) XMERIT*Oa 
Дино у тобт еј еј) XMERITS1v 


melt (ls ia NELEM)EXMERIT 


ХЧЕР | Т + ( |" ( ХНСМ (Т.)) SFMEANSHO0 05) D¥( Lew TMee(leeTM1) ) 


ЕТ. Т2.) XNERIT =0 
12 (ХМЕВ1Т.ОТ.1.) XMERTT=le 


ES Gear te NELEM)=xXMER] T 
ИСЕ ХАРЕ О: 

ПЕРУ СЕ»9) G8 T8 301 
RETURN 


E SN Е ОАК ТТУ РОК USE 
шас ыва 200 TIS USEC FOR 


O OTIS A 


вен 501 Теје ем 
БО: 1351,2 
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INFINITY 


Settee OE CODER 





154: DESSERT Toy ESO) MERTT(11,13,12)2200+; 


1:551 haeo 501 

155; ШК 3 [2 -е-АГ СОСПЕКТІТ(І1,13,12)) 
MS OL CCE Е 

132: EIER ее Еі 

С 

160: Г и ЕЕ 1900$ FOR NEXT ITERATION 
Je l C 

ше зи 08 600 1-1,11 

164: МАМЕА ОР ЈОРАСЕ (1) 20 

Moon EOC CONTINUE 

166: БО» 11=[/4 

167; Some do 12-115 

ToS; 072.05 13=132 

EO CG TITTIES 

I NELEM=0 

17 1: RETURN 

шие: ЕМО 
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ЕКО BF АРН БОТ E 


p> 
О «О WNC gl 


— pa ja 
ој го => 


14: 
=; 
or 
7: 
фа: 
i 
С; 
ea? 
га, 
23: 
25: 
Bo 
26: 
27: 
B 
B 
mac 
31: 
Bo: 
5 E 
341 
25: 
36: 
37: 
38: 
85: 
4O0: 
41; 
42: 
43; 
44: 
45: 
46; 
à 


BOSE SUTINECDECODER 


flow overcomes TePLEMrENTS THE 
ОКР 


C; Cy Gy cy G 


REAL LSAV,MERIT,LAMDA,LENSTH, NU 

REAL LNSTHI>LNSTH2ZSILNGTH3,LNGSTHU, LNGTH5, LNGT 46 
INTEGER SEQSAV,CHAR 

КОО AUC Те ), 1оаме Ссбе»12), 5 АМЕ (5) 

ПО и (11), ХТКА5 (11) 

DIMENSTI@N LENGTH(6),YLNGTH(6)s2nLNGTH(6).,ZLNGTH(6) 
DIMENSISN MERIT(4,2»12)»CHAR(6512) 

DNE ONSNAEMATKUIT)A,NHSPACEC1II1) 
DIMENSDISNONSEPCUCT2)ANSQ(12)^SEGSAV(S5, 12) 

АЕ ДАСКЕ МЕ ЕМ, МЕКТТ 

С Г УВИЕСКЗ/ТРАМО 
СОММЗК/В| еСК4 / МАК, КОРАСЕ, МСНАК 

ПАЛА Л ОУ 


INITIALIZATION 


OOO 


О8 50 1517128 
50 NSP OS ЈСО 


Uc EUM SDI CITIES GF EACH CHARACTER 
КЕКЕ О ОТА ЛАТ ОЧ 


СҮ СУУГУ 


LENGTH( 1) s"AL8G(0+269) 
LENGTH(2)s-AL8G(O 351) 
ШЕ -ІЕРСТЯ( Уа-4Д| 95(0.159) 
eign (>) =LONCTH(G)=-AL9G(0«232) 
ОВО Т 120 
МАТА )=ZENGTH(1)=YUNGTH(])=e00- 
120 CPNTIRUE 
Селище 151,11 
XTRANS( J) 206 
TOT СЕТЕ 
ОЗИ [126 
рама ск а 12 
LSAV(I4J)20* 
[ЗАМЕ (] 14) ЕО 
10 CoS TINUE 
С 
С Ae Ss USED TO SELECT OUTPUT OPTION 
С 
IF (SENSE SWITCH 2) 8,9 
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4R: 
49; 
SOS 
917 
See 
s 
54: 
ser 
56. 
57: 
23: 
ба. 


= = Y 


60% 
6 1c: 
ба, 
аз: 
64: 
бЭ 
66: 
OU 
68: 
6. 
79: 
До 
деи 
8: 
74: 


FOTO TED 


Е 


@ 
С 


5 WRITE(6,3C00) 

3000 FORMAT('1") 
ARITE (6,4010) 

4O10 FBRVAT(SX,'!SURVIVS8R SEQUENCES!) 
Neem OC) (ly l=t,il) 

ВЕ te (1 cs 5X ) ) 

> CONTINUE 


MAIN ALGERITHM 


Е: 15 А ОТ 

ПЛО те ANCSEEMCNT SPACE 
Е Е LETTERS PACES 
Е Ө ARE DASHES 


NELMI=NFLEM+] 

DO 90C KalsNELMI 
Пе ы) GE TS S10 
CO TORTE] 


ODO SIT ISIS] 


XTRANS (J) STRANS (CY) 
ERU МТ: Е 
В Е» GO Te 920 
possc sats 14 
XTRANS (J) 208 


СЕТІП ОЕ 


ПОВЕ І 

LAMDAZLENGTH(2)*XTRANS(3) 3 Jee 
NUSLENGTH(3)4XTRANS (5) 
т AIDA) САМБА ЧОТ ]1=3 
NUFZLNGTH(4)4XTRANS (5) 
IF(NUSLT«LAMDA) LAMDAzNU 3 [54 
ISAVE (12k) =] 
ENGTHISLAMDA+MERITT( isla) 

300 CONTINUE 


NSDE 2 

LAMDA=LENSTH(1)+XTRANS(1) 3 ]=1 
NU= ZLNGTH(%+)+XTRANS(6) 
IFINUYeLTeLAMDA) LAMDASNU 3 l=4 
NU=LENGTH(5)+XTRANS (9) 
IM LT LADA) LAMDA=NU 1 185 
NU=SZLNGTH(6) +XTRANS (9) 
IF (NUeLTeLAMDA) LAMDA=NU ; 126 
ISAVE (Ps kK) =] 
LNGTH2sLAMDA«MERIT(2, 15X) 

Son NE 


МӘСЕ S 


Pete men oly ХТ Ал 512) ) 141 
NUSLENGTH(5)+XTRANS (10) 
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ques. 
10897 
199% 
104; 
тави 
106: 
10% 
1027 
109% 
LEOS 
1 5 
КЕЕ 
По 
114: 
mls: 
Ls 
117 
та 
1192; 
120. 
m 
lo 
low 
Tea 
аб 
Шабе 
er. 
Leo. 


1308 
т. 
io 
13% 
e 


155: 


136: 
jo 
ju 
139: 
140: 
14 
1421: 
143: 
144: 
145: 
146: 
1473 
148: 
149; 
ТЕЛІ; 
151. 
qo 
Ша зе 


62 
С 
E 


304 
С 
С 


BCS 


С 
С 
С 


14 


ИИ Е ЛЕВ отне FOR SURVIVOR 


ІБІК РІ Т»ҺАМДА) ГАМОАЕКИ ; [25 
NUSZLNGTH(6)4XTRANS (10) 
рамо теб АЧРАЗ C EAMOATNU 3 
NUZLENGTH(Z)+XTRANS(11) 

A ОСА) LAMOASNU 3} 
ISAVE(3,X)=1 

"МОТНЗ«| АМОАЗМЕЗІТ(3,1,%) 
CONTINUE 


I =° 


AOE А 


MANDA ER TEA (IE XTRANS(2) J ТЕЦ 
NUSZ7LNGTH(4)4X TRANS (7) 

ЈЕ (МИ• Те АМА) LAMDAZNU ; 
NUZLENGTH(5)4XTRANS (14190) 

J (U S Т.БАМОА) ГАМВА=МИ 3 
БОЕ E GTH (G)+X TRANS (10) 
Maeve sl AOA) БАЗБАВ МИ 3 
NUSLENGTH(2)4XTRANS (11) 
ТЕТ. ВАА) LAMSASNU 3 
ISAVE (4sK) =] 
LNGTH4SLAMDAFMERIT(3,20%) 
CANTINUE 


ШЕШЕ, > 


ESL -PENGTHI2)*XTRANS(A) ; 
МО=СЕКСТН (3)+ХТКАМЗ (8) 
ша иза МОЛУ ГАМОА=МИ 3 ]=3 
NU=SZLNGTH(4)+XTRANS(8) 

ТЕ АМИ ЕТ. САМА) LAMDAENU 3 
ISAVE (5,«%) 2] 
LNGTHSszLAMDA«"ERIT(A4, 12K) 


CONTINUE 
М8СЕ 6 
LAMSDA=LENGTH(2)+XTRANS(4) ; ]=6 
NUSLENGTH(3)4XTRANS (8) 
1Е (МИ. СТ. САМОА) АМГА = МО $ 1+3 


NU=ZLNGTH(4)4XTRANS (8) 
JF(NULT.LAMDA) LAMDAZNU 3 I=% 
ISAVE (6,K)=1 
LNGTHO6sLAVDA*ME2IT(4,2.K) 
CONTINUE 


SE амезе 


BE cO Si ЕНІ 
LENGTH(?2 ysUNGTH2 
CENA TREO =LNIGTAS 
LENGTH (4) =LNGTH4 
аза аа а Це ТН 
LENSTH(6) =LNGTHE 

пр па 41,6 

Bora CENGTH(I ) 
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186! 
ray: 
Шери 
1а 
rgo: 
po 
1225 
193% 
194; 
eo oe 
196: 
197 
198: 
199% 
200% 
coe 
о 
Coe 
eov 
205“ 
COG: 


Са Сэ оез 


(> Cy C OO 


200 
EO 


DEN200 ow Ey 6 

Ze а Мо ТАСТ); 
ПА руга матн (1) 
DUINCOINHCIOEUENGTHLUUDT) 
СВОЕ 

CONTINUE 


DETERMINE MINIMUM LENGTH SEQUENCE FROM THE 
STRESURVIT VSR SERUENCES 


2091 


са 


2i 


So 


40 


ГА) 


73 


NUFLENGTHCL) 5 
DO 801 I=2»6 
Е pU 
CONTIAYE 

WE IGHT=LENGTH(ISURV) 
1Ме | 

м= 0 

NETO AVE (TOSURVI NELI) 
ПРЕСОТ ЕМ 

CONTINUE 

IFCON SEQ 4) efRe (No EC eG) ) 
МЕМ+ | 
ТЕССКЕГМІ-М).(Е.»90) 
NS= ISAVE (Na NELM1=M) 
Gu To 31 

МЕеМ З 

ПР СЕ ем) е | Е 40) 
NOSISAVE (NA ANELMI= M) 
ЛЕН ў 
NSEQ(C IM) =&NENS 

Se ја ве 

BO NTIBUE 

09 40 151,1М 
NSQ(1)=NSEQ(IM-1+1) 
L=L+1 
NSAVE(L)ZzWEIGHT 

ОО Л psu pM 
SEQSAVCIEXI)*NSG(CI) 
1М 1 + | М + 1 

DMA ТИТхІе 
SEQSAV(L, 1) #0 


ISURV= 1 


ЕТ Я ISURV=] 


an) 


CU 9.35 


се 19 


Ош РОТ 


sum Ul OPTION 1 


60 


T2 
4001 


MENSE SWITCR 2) C0270 


weet LEO УТРУ? 
CONTINUE 

Е гб 

КИТЕ ба 001) |), (ТЗБАУЕ(.,<)» Кс 19 ЧЕ( ЕМ) 
аа Le (11+6X)) 

АРІТЕ(6,4002) 


146 





ЕШ: 
сос; 
209: 
со: 
с: 
Ера: 
EIS 
214: 
се: 
216: 
217: 
218; 
219: 
соо: 
се! 
вес: 
RRS 
ру: 
ВЕР: 
‚се: 
Ber, 
een. 
со; 
E30, 
al; 
езг: 
es), 
C34; 
AO y 
ЕС 
eo 
238% 
Е 
240; 
241: С 
све. С 
242% С 
244} 
245, 
246: 
CAT: 
pus 
ЕН 
ага е 
== 
сы: 


wii. Ф 


258; 


(O02 РВКЕМАТСИ) 

ма 1ТЕ (6,5000) 
SOO ИКО SURVIVOR LENGTHS?!) 

раз 221%6 
13 WRITE (624003) Je (LSAV(JUsK),K=1,NELEM) 
ПІЗ БӘӨПКУЕТТЕХУ, | Пусха 12 (бег, 1 Х)) 

WRITE (624002) 

WRITE (621003) 
DOCS pane meee ET KEL THSeCS! ) 

1.2221 

Ге 

ВА 15992 

Jar, 
Dee) Jly(3ERIT(1,lJ iX) I K=1, ,NELEM) 

08 9205 Is3,4 

Бе) 71.2172 

А а 
ЛТ Јо С ПЕРЕТСОЈРТЈК) Кај» МЕЦЕМ) 
ОСА Ра а Ро <= ОТ ах ЦС (Ебесол 1 Х)) 

WRITE (56/4002) 

WRITE (6,4021) 
ASS Xx) DECIDED?) 

ARITEC6+ 4020) 
ОСО ха МАЕК ЗХ, "ӘРАСЕ!, 1 COX, "ПІ ЕКОТЯ!,9Х, 

II SEEN?) 
90% wWwRITE(621000) MARK(1),NSPACE(1)-WEIGHTs (NSQ(K), 

кем) 
TOCO ШО Ха | > 3х, | 2» ] 1Хх»Р794 3,107», 1211) 

D3 96C Is2,NCHAR 

AMENS 1000) MARK(I)NSPACEt(I) 
SEO CONTINUE 

ЕЦ јез) (20 

свары 72 
726 ие.) 69 та 72 


Oe eee ea oN eee DECODED SEQUENCES AND LENGTHS ONLY 


(29 
БЕТШТЕ(6,6000) (АБАМЕ(1]1)»121»5) 
s I (66001) ISE ISAV (TJ) Ss 2012), 1=21,0) 
с о х го ОЕ 9. 3л2х) 57) 
COUT аза 1 о х 20 (1111»1%Х)) 
WRITE (626003) 
6003 FORMAT (//) 
үс. СӘМТІММЕ 
KETURN 
END 
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AFORTRAN 
dg: 
с 
er 
4; 
E 


NM ON оо Y 0 YU + Y (Q > O xD бо ~ 


41: 
42: 
43: 
44: 
45: 
46: 
47: 
481 
491 
5б 
SL. 


LSsG8 

АЈА ВЕЦ THE PRE-CETECTION FILTER 
C DESCRIBED IN THE THESIS TEXT 

ё 


REAU CDOT MEANX; IDENT 
IESO TON ре Т Сг, е), АІ а), н(е) 
BETEN ENTE г), ТЕГРЧ1 (2.2), ТЕМРЯА (2), ТЕЧРС (2) 
ОИУ ИО ЕВ 23200 РЧАВ (2,0) РНІ (2,2) х РНІТ (2,02) 
DIMENSION ISUF (4000) 
Мат От LOST, OUASH,T»TI¢gTesFREG 
ӘРТІСІ) "ENTER T AND FREQ*! 
ОО) 'LDET AND DDASH' 
10) 
E 
E INITIALIZE VARIABLES AND INDICES 
С 
ХР1=ХР2=0.5 
ТАЦ=0.00025 
KT=0 
KIND=JIND=0 
XBUTsOs 
МЕАМХ=О. 
5520. 
МАВХ= 0.5 


Тэд. Т i LDO8Ts^«*LOGT; DDASHz4s*DDASH 
11=4.*71 ; То=4. «Т? 

ҒКЕ?г6»28319ҡх<ЕҒрЕ0 

1023«•+1 

тво 


20 CONTINUE 
PH1(112)=PH1T(2+1)=FREG*TAU 
ӘБРІССО, 1) <РИТТ(1,2):-ЕрЕФ«ТА)() 
Du 151/2 
Пон Jete 
Пра ТО ја југо. 

EVAR(IG JJ)s0* 
(Тл) =0 

el CONTINUE 
АДА РА У ХРМАК( 272) =1• 
ПОТ >» 1 )=1рЕ ТТ (272) =1• 


Н(1):1. 
н(23)-С0.» 
100 X309 NT TNIUE 
8 
С INPUT DATA FRB“ TAPE 
С 
EB c ei cUEPERTN (1) 1s 1 RUF 24000, TERR) 
1 Ди вв Е Це 1) СЕ 19 1 
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за: аеш аша, Зад Да) ТЕКЕ 
СОСТ S O99 UINDEXsTI,*4000 


54: С 

85% С BEGIN PROCESSING 

55C 

S^. e zuEBATCOUSUPCINDEX))/2*223 
55: Хеви] 

59: Х1=РиХ] 

60: ЈА O 1000009) KT=O 

с: КТЕКТ + | 

en CACIT STATS 

oc. С 

64: C GAIN 

65: CALL MVMULTIPVA?R ¿md TEMPC) 
66: CALL VVMULT(M,TEMPC, TEMP) 
2728 | ТЕМО 15 |е/ ( ТЕМР + МАРХ ) 

os: CALL MVMUILTI(PVAR,H,TEMPC) 

сз И ЕТЕ РС (] & ТЕНР] 

p СА1М (2) =ТЕМРС (2) „ТЕмР] 

Peo ie 

p ESTIMATION 

из: XH1=xP1+GAIN(1)x (X1-xP1) 
74; XH2=XxXP2+GAIN(2)* (X1-XxP]) 
p С 

SNC Шаша Ка СЕМ VARIANCE 

3 TE mop У ИДО ТОЗАЛ МУНА ТЕМРМ) 
78: ВОО 1512 

EE рис ОС ШЕ Та Е 

ар ЕИО = Трем ТА у ТЕКРМ (ІА) 

81: ШӘШЕ ГТГТЕ“РМ;:РУАПН,)ГУАхя) 

С 

E PREDICTISN VARIANCE 

84: CALL MMULT(EVARsPHIT,TEMPM) 

EE СА ПИПІГГІРЫІ»ТЕМРЖ,РУАК) 
во С 

СС ВН ЬТЕР PREDICTISN 

88: ХР | 5ХНЪ+РНТС | )Р) «ХЧ2 

So. ХРе=РрРН]1 (2,1) «ХН1 +ХН2 

90: UC 

С SeUSke FPILTERES ESTIMATE AND LOWPASS FILTER 
eC 

ENT Х=ХН1**жё 

Eu XK T =K T 

95: Ме ЮТ»С5Т50) G69 T9 60 

96% ASUS XSGUT+ (1. /XKT)*(X-X2UT) 

E G9 T8 61 


cree OO XBUTEXGUT+O «OP? * (Xe XGUT) 
2.406) CANTINVE 


КОС; [ХНАТ=-] 

IO ПАСУ ез е МЕ ЛМХ) ПХНАТ Е 1 
102; CALL МАЗИ 

TO е 


104: С СЕИ APUESTO D/A ROUTINE FOR ANALAG RECORDING 
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105% 


AND 


(СЗ СУ СУ СУ 


ХХНАТЕР | ВАТ (ТХНАТ) 
CALU DAL((X2^X15XH15 X8UT^ GAINCLIZEVAR( Ts 1s EVAR (Cre), 


1МАКХ,Р,ХХНАТ) 


CONTINUE 
Goya Tee 


РЕ ИЕ STATS 


БЕТЕР ESTIMATES THE NOISE VARAINCE 


XKT=KT 

IF (КТ+СТ, 4000) G8 TE 70 

МЕАМХ= МЕАМХ+ (1. ИХКТ ) # (ХОЏТУ-МЕАМХ ) 
Gans 7] 

МЕАМХ ЕМЕАМХ + (16 74000 • ) У (ХВЏТУМЕАМХ ) 
CONTINUE 

IAS. GT. “EANX) G8 T9 80 
SoeootoeQix(XGUT=SS) 

БӘЛІШІ ШЕ 

VARX=SS 
IF(VARX«LE:O.O]J) 
RETURN 

END 


VARX=0.O1 
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АҒӨЗТКАМ (5,69 
ie, SUBRBUTINE PMULTCA,B,C) 
С Miler two MATRICES 
(И НМ А се гејав се» 2)»С0 02» 2) 
Шао тетла 
ОВ 10 Ј=1,2 
muero) so. 
РАНО Те рус 
D3 20 J=1,@ 
OA Уг 
СЕТА Ј) =А(1К) +8 (Ки .Ј)+С(1/ 7) 
CONTINUE 
RETURN 
END 


هم 
O‏ 


Y 
с? 


е ве ее oo зо е е "o e фо ве зо зе 


u) Г + CO «0 00 - € Ul ur 


هم هم همهم 
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END OF COMPILATION 

1: SUSE ИДО Т СА, ВС) 
С ПВ TR IX BY COILUMN VECTOR 
Ole STEN ACerc)sd(e)sCle) 
ро Fe Mi =) 2 
ЕСТЕ Ы 
C9 20 1:10 
DER Іа 
CCI)sACI,K0*B(K)4CÓI) 
CONTINUE 
RETURN 
END 


ph 
O 


e C) Оо NO УО Ul г шг) 
го 
о 


ээ ес ее эъ зе оо зе ае vo. Фе 


هم هم 
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ANI 


هم حم 


r o pil AT van 

SUBRSUTINE VMMULT(AS, B.C) 
C Ищар ОО ЧЕСТОВ BY MATRIX 
ВЕ а Аса да (деж, 2), Се) 
DO IO ISIC 
зе Со 
ео 151902 
шаша Ке | хуЕ 
БОРА КВ (К 1)+С 01) 
CONTINUE 
RETURN 
END 


> 
O 


а м + ہہ (۸ (ر)‎ 
го 
= 


оо ге ев ее ге аа Фе ee о ее езе 
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BP) OF PCH Riad DEN 
1: SUSREUTINE УУМОГТСХА,В,С) 
С Ем УСТЬЮ BY COLUMN VECTOR 
ОЗ БОМИ Cela Sle) 
СЕО, 
ОЛО Tele 
СЗА (|) В (Т) +С 
10 CONTINUE 
RE TURN 
END 


VO CONE DNS WM 








END 8F COMPILATION 
1 РОТЕ м (АБ. С) 
С MULTIPLY CSLUMN VECT3R BY ROW VECTOR 
OTE EISSN А(27,8(27,СС(2,2) 
МАО Теје 
(ОКО 2751,2 
АТВ) 
CONTINUE 
RETURN 
END 


мо OO Y DU + uy ГО 
pa 
O 
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